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Abstract 

The  development  of  a  quantum  computer  based  on  a  system  of  trapped  atomic  ions  is  described, 
following  the  proposal  of  Cirac  and  Zoller.  Initial  results  on  a  two-bit  quantum  logic  gate  are  pres¬ 
ented,  and  select  experimental  issues  in  scaling  the  system  to  larger  numbers  of  ions  and  gates  are 
treated. 
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I.  Introduction 

In  its  simplest  form,  a  quantum  computer  is  a  collection  of  N  two-level  quantum  systems 
(quantum  bits)  which  can  be  prepared  in  an  arbitrary  entangled  quantum  state  spanning  all 
2n  basis  states  [1,  2],  A  quantum  computer,  unlike  its  classical  counterpart,  can  thus  store 
and  simultaneously  process  superpositions  of  numbers.  Once  a  measurement  is  performed 
on  the  quantum  computer,  the  superposition  collapses  to  a  single  number,  which  in  some 
cases  can  jointly  depend  on  all  of  the  numbers  previously  stored.  This  gives  the  potential 
for  massive  parallelism  in  particular  algorithms  [3],  most  notably  an  algorithm  which  fac¬ 
torizes  numbers  efficiently  [2,  4],  Apart  from  applications  to  this  and  other  algorithms  [5], 
creating  multi-particle  entangled  states  is  of  great  interest  in  its  own  right,  from  the  stand¬ 
point  of  quantum  measurements  [6]  and,  for  example,  for  improved  signal-to-noise  ratio  in 
spectroscopy  [7,  8]. 

Unfortunately,  there  are  very  few  physical  systems  which  are  amenable  to  the  task  of 
quantum  computation.  This  is  because  a  quantum  computer  must  (i)  interact  very  weakly 
with  the  environment  to  preserve  coherence  of  the  superpositions,  and  (ii)  interact  very 
strongly  with  other  quantum  bits  to  facilitate  the  construction  of  quantum  logic  gates  neces¬ 
sary  for  computing.  In  addition  to  these  seemingly  conflicting  requirements,  the  quantum 
bits  must  be  able  to  be  controlled  and  manipulated  in  a  coherent  fashion  and  be  read  out 
with  high  efficiency. 

In  1995,  Cirac  and  Zoller  showed  that  a  collection  of  trapped  and  cooled  atomic  ions 
can  satisfy  these  requirements  and  form  an  attractive  quantum  computer  architecture  [9].  In 
their  proposal,  each  quantum  bit  is  derived  from  a  pair  of  internal  energy  levels  of  an 
individual  atomic  ion.  By  using  laser  beams,  the  quantum  bits  are  coupled  to  one  another 
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by  virtue  of  the  quantized  collective  motion  of  the  ions  in  the  trap,  mediated  by  the  Cou¬ 
lomb  interaction.  A  reduced  version  of  their  scheme  was  implemented  in  an  experiment  on 
a  single  trapped  ion  [10]. 

This  paper  concentrates  on  the  trapped  ion  quantum  computer  architecture,  and 
covers  some  of  the  experimental  details  involved  in  conducting  simple  logic  opera¬ 
tions  between  small  numbers  of  trapped  ions.  Following  a  brief  introduction  to  ion 
traps  and  the  interaction  between  internal  and  motional  states  of  trapped  ions  in 
section  II,  preliminary  experiments  involving  one  ion  are  reviewed  in  section  III,  and 
particular  technical  problems  in  the  extension  of  this  scheme  to  N  >  1  ions  are  dis¬ 
cussed.  Section  IV  considers  the  problems  of  addressing  individual  ions  in  a  string 
with  lasers  and  the  cross-coupling  of  the  3 N  quantized  modes  of  motion,  and  section 
V  characterizes  some  expected  sources  of  decoherence.  The  topics  covered  here  are 
by  no  means  exhaustive,  but  may  give  an  indication  of  some  of  the  key  problems 
which  may  lie  ahead  in  the  near  future.  We  have  attempted  a  more  complete  investi¬ 
gation  of  the  problems  in  Ref.  [11].  For  more  reviews  of  the  trapped  ions  quantum 
computer,  see  Refs.  [12,  13].  Other  notable  physical  systems  proposed  for  quantum 
computation  not  covered  here  include  cavity-QED  [14]  and  bulk  spin  nuclear  mag¬ 
netic  resonance  [15]. 


II.  Background 

A.  Internal  states  and  detection 

Ions  can  be  confined  for  days  in  an  ultra-high-vacuum  environment  with  minimal  pertur¬ 
bations  to  their  internal  atomic  structure,  and  collisions  with  background  gas  can  be  ne¬ 
glected.  Even  though  the  ions  interact  strongly  through  their  mutual  Coulomb  interaction, 
the  fact  that  the  ions  are  localized  necessarily  means  that  the  time-averaged  value  of  the 
electric  field  they  experience  is  zero;  therefore  electric  field  perturbations  are  small.  Mag¬ 
netic  field  perturbations  to  internal  structure  are  important;  however,  the  coherence  time 
for  superposition  states  of  two  internal  levels  can  be  very  long  by  operating  at  fields 
where  the  energy  separation  between  levels  is  at  an  extremum  with  respect  to  field.  For 
example,  a  coherence  time  exceeding  10  minutes  between  a  pair  of  9Be+  ground  state 
hyperfine  levels  has  been  observed  [16].  It  is  also  possible  to  employ  a  ground  and 
excited  (metastable)  electronic  state  of  a  trapped  ion  as  a  quantum  bit  [9].  This  option 
seems  difficult  at  the  present  time,  primarily  because  the  energy  splitting  is  typically  in 
the  optical  region,  requiring  extremely  high  laser  frequency  stability  to  drive  coherent 
transitions. 

Figure  la  shows  a  reduced  energy  level  diagram  of  a  single  9Be+  ion.  Although 
many  other  ion  species  would  also  be  suitable  for  quantum  computation,  we  will 
concentrate  on  9Be+  here  for  concreteness  and  to  make  a  connection  to  the  experi¬ 
ments  at  NIST  [10,  17].  We  will  primarily  be  interested  in  two  electronic  states,  the 
2Sm(F  =  2,  /(if  =  2)  and  2Si/2(F  =  1,  mp  =  1)  hyperfine  ground  states  (denoted  by 
||),  and  ||),  respectively),  separated  in  energy  by  Hcoq.  These  long-lived  “spin”  states 
will  form  the  basis  for  a  quantum  bit.  Detection  of  the  spin  states  can  be  accom¬ 
plished  using  the  technique  of  quantum  jumps  [18].  By  tuning  a  polarized  laser 
beam  to  the  |X)  — 2 P3 /2  transition  near  313  nm  (Fig.  la),  many  photons  are  scat¬ 
tered  if  the  atom  is  in  the  ||)  spin  state  (a  “cycling”  transition),  but  essentially  no 
photons  are  scattered  if  the  atom  is  in  the  ||)  spin  state.  If  a  modest  number  of 
these  photons  are  detected,  the  efficiency  of  our  ability  to  discriminate  between  these 
two  states  approaches  100%. 
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(a)  (b) 

Fig  1 :  (a)  Electronic  (internal)  energy  levels  (not  to  scale)  of  a  9Be+  ion.  The 
2S[/2 (F  =  2,  nif  =  2)  and  2Si/2(F  =  1,  nip  =  1)  hyperfine  ground  states  (denoted  by  ||) 
and  |t)  respectively),  separated  in  frequency  by  (Oq/Ijt  m  1.250  GHz,  form  the  basis  of  a 
quantum  bit.  Detection  of  the  internal  state  is  accomplished  by  illuminating  the  ion  with  a 
cr+-polarized  “detection”  beam  near  313  nm,  which  drives  the  cycling 
2S1/2(F  =  2,  Ttif  =  2)  — >  2P3/2 (F  =  3,  nif  =  3)  transition,  and  observing  the  scattered 
fluorescence.  The  excited  P  state  has  radiative  linewidth  y/2jr  a  19.4  MHz.  (b)  Energy 
levels  of  a  trapped  9Be+  ion,  including  the  motional  states  of  a  single  mode  of  harmonic 
motion,  depicted  by  ladders  of  vibrational  states  separated  in  frequency  by  the  mode  fre¬ 
quency  a>z.  Two  “Raman”  beams,  both  detuned  A  to0>  o)z,  y  (  from  the  excited  2Pl/2 
state,  provide  a  coherent  two-photon  coupling  between  states  | n)  |J.)  and  | n')  ||)  by  setting 
the  difference  frequency  u>l  to  match  the  desired  transition  frequency.  As  shown,  the  Raman 
beams  are  tuned  to  the  first  red  sideband  of  the  ||)  — >  |T)  transition  (qjl  =  co0  —  a>z  ). 

B.  Ion  traps  and  motional  states 


In  Fig.  2,  we  show  a  schematic  diagram  of  a  linear  Paul  trap  [19],  consisting  of  four  elec¬ 
trode  rods.  The  linear  trap  is  similar  to  a  quadrupole  mass  filter  [20]  which  is  plugged  at 
the  ends  with  static  electric  potentials.  A  potential  Vo  cos  is  applied  between  diagonally 
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Fig.  2:  The  upper  part  of  the  fig¬ 
ure  shows  a  schematic  diagram  of 
the  electrode  configuration  for  a 
linear  Paul-rf  trap  (rod  spacing 
~1  mm).  The  lower  part  of  the 
figure  shows  an  image  of  a  string 
of  199Hg  ions,  illuminated  with 
194  nm  radiation,  taken  with  a 
UV-sensitive,  photon  counting 
imaging  tube  [22].  The  spacing 
between  adjacent  ions  is  approxi¬ 
mately  10  pm.  The  “gaps”  in  the 
string  are  occupied  by  impurity 
ions,  most  likely  other  isotopes  of 
Hg+,  which  do  not  fluoresce  be¬ 
cause  the  frequencies  of  their  re¬ 
sonant  transitions  do  not  coincide 
with  those  of  the  194  nm 
2Si/9  P1/2  transition  of 

1 QQ,'.  i  ' 


opposite  rods,  which  are  fixed  in  a  quadrupolar  configuration,  as  indicated  in  Fig.  2.  To 
provide  confinement  along  the  axial  z  direction,  static  potentials  Uq  are  applied  to  the  end 
segments  of  the  rods  as  indicated.  We  assume  that  the  rod  segments  along  the  z  direction 
are  coupled  together  with  capacitors  (not  shown)  so  that  the  rf  potential  is  constant  as  a 
function  of  z.  Near  the  axis  of  the  trap  this  creates  a  potential  of  the  form 


0 


Vo 

2 


cos  Qjt 


U0  (2 z2-x2-y2\ 

2  V  d2  )’ 


(1) 


where  R  is  equal  to  the  distance  from  the  axis  to  the  surface  of  the  electrode  and  d  is  a 
characteristic  axial  dimension  of  the  static  electrodes.  This  gives  rise  to  harmonic  pondero- 
motive  potentials  [21]  in  the  radial  ( x  and  y)  directions  accompanied  by  static  harmonic 
confinement  in  the  axial  (z)  direction  resulting  in  an  effective  3D  harmonic  confining  po¬ 
tential  U(r)  =  |  m corx2  +  i  mcj2y2  +  i  mcu2z2,  with 


(2) 


where  cop  =  qVo/ (21'2  QjinR2)  describes  the  ponderomotive  portion  of  the  potential,  and  q 
and  m  are  the  charge  and  mass  of  the  ion,  respectively.  In  these  expressions,  it  is  assumed 
that  mp  <  Qt,  a  condition  known  as  the  “pseudopotential  approximation”  [21]. 

Figure  2  also  shows  an  image  of  a  “string”  of  I99Hg+  ions  which  are  confined  near  the  z 
axis  of  the  trap  described  in  Ref.  [22].  This  was  achieved  by  making  cox,  coy  Zz>  c»z.  thereby 
forcing  the  ions  to  the  axis  of  the  trap.  The  spacings  between  individual  ions  in  this  string 
are  governed  by  a  balance  of  the  force  along  the  z  direction  and  the  mutual  Coulomb 
repulsion  of  the  ions.  Example  parameters  are  given  in  the  figure  caption. 

Of  the  3N  normal  modes  of  small  oscillation  in  a  linear  trap,  we  are  primarily  inter¬ 
ested  in  the  N  modes  associated  with  axial  motion.  A  remarkable  feature  of  the  linear 
ion  trap  is  that  the  axial  modes  frequencies  are  nearly  independent  of  N,  offering  the 
possibility  that  mode  interference  might  be  small,  even  for  large  numbers  of  ions  [9].  For 
two  ions,  the  axial  normal  mode  frequencies  are  at  ai7  and  \/3o)-;  for  three  ions  they  are 
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( oz ,  \/3  a>z  and  \/5.4  a>z.  For  N  >  3  ions,  the  Mh  axial  normal  mode  can  be  determined 
numerically  [12,  13,  23].  We  will  concentrate  on  the  axial  center-of-mass  (COM)  mode, 
which  is  the  lowest  frequency  mode  and  is  also  the  most  resolved  in  frequency  from  the 
others.  The  quantum  state  of  axial  COM  motion  at  frequency  a>z  can  be  described  by  the 
ladder  of  vibrational  eigenstates  | n),  of  energy  1icoz(n  +  \)  with  vibrational  index  n  de¬ 
scribing  the  number  of  “phonons”  contained  in  the  collective  harmonic  COM  motion. 


C.  Coupling  between  internal  and  motional  states 


We  describe  the  coupling  between  the  internal  levels  of  a  particular  ion  in  the  string  with 
the  COM  axial  mode  of  collective  motion  when  a  classical  radiation  field  is  applied  to  that 
ion.  If  the  internal  levels  |J.)  .  and  |f).  of  the  yth  ion  in  a  string  are  coupled  by  a  dipole 
moment  operator  / *■  (other  couplings  can  be  shown  to  behave  analogously),  then  exposing 
this  ion  to  a  traveling-wave  electric  field  E(r)  =  Eq  cos  (k  ■  r  —  a>it  +  qp)  with  frequency 
a>L,  phase  cp,  and  wavevector  k,  results  in  the  interaction  Hamiltonian 

=  -fij  ■  E[r )  =  h£2>\ [Sj+  +  Sj_)  (e'Wz-^'+'P)  +  e-i{k  z-<oLt+v ))  _  (3) 

In  this  expression,  Q1  =  —  (T|  fij  ||)  -Eq/ATi  is  the  resonant  Rabi  frequency  connecting  |J,)  . 
to  |f),  in  the  absence  of  confinement,  S'+  (S'_ )  is  the  internal  level  raising  (lowering)  opera¬ 
tor  of  the  jth  ion  changing  ||)  •  to  |f)  ■  and  vice-versa,  z  =  Zo z{ae~”°z<  +  a' ew>zt)  is  the  axial 
COM  coordinate  operator  of  the  confined  motion  with  associated  harmonic  raising  (low¬ 
ering)  operator  a)(a )  and  zero-point  spread  zo  =  (h/2Nm(oz)] ,2 ,  and  Nm  is  the  total  mass 
of  the  ion  collection.  If  the  applied  radiation  frequency  is  tuned  to  <x>l  =  a> o  +  (n1  —  n)  <»z, 
thereby  coupling  the  states  | n)  ||).  and  | n1)  |f)-,  Jf7-)  is  transformed  to 

=  h£V{$+  +  S>_  (4) 

in  a  frame  rotating  at  mL,  where  terms  oscillating  faster  than  Q1  (Q1  w o)  have  been 

neglected.  Here,  rj  =  (k  ■  z)  "n  is  the  Lamb-Dicke  parameter,  which  controls  the  amount  of 
coupling  between  internal  and  COM  motional  states. 

When  the  coupling  of  Eq.  (4)  is  applied  between  the  /th  ion  and  the  COM  mode  of 
motion,  the  system  evolves  between  the  two  quantum  states  | n)  |f).  and  |n')  |f).  with  Rabi 
frequency  [24,  25] 


=  O’ \(ri\  e'7(fl+flt)  |„)|  =  Qijj K-»l  e^1!2  x  , 


«>  ! 


(5) 


where  «>(«<)  is  the  greater  (lesser)  of  n  and  n'  and  (x)  is  a  generalized  Laguerre  poly¬ 
nomial.  If  the  Lamb-Dicke  criterion  is  satisfied,  where  the  amplitude  of  the  ion’s  motion  in 
the  direction  of  the  radiation  is  much  less  than  1/2 jt  (or  nl/2t]  <  1),  we  can  evaluate  Q’n, 
to  lowest  order  in  17  to  obtain 


Q> 


\n'-n\ 


—  n 


(6) 
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We  will  be  primarily  interested  in  three  types  of  transitions  —  the  carrier  («'  =  n),  the  first 
red  sideband  («'  =  n  —  1),  and  the  first  blue  sideband  (, n '  =  n +  1)  whose  Rabi  frequencies, 
in  the  Lamb-Dicke  limit,  are  given  from  Eq.  (6)  by  O’,  and  rj(n  +  1 ) 1  /  2  Q1  respec¬ 

tively. 

In  practice,  driving  direct  transitions  between  | n)  ||).  and  j «')  ||).  with  rf  or  microwave 
radiation  is  not  be  feasible,  as  the  sideband  operation  transition  rates  (proportional  to 

=  2nzo /A)  would  be  extremely  slow  due  to  the  long  wavelength  of  the  radiation.  Alterna¬ 
tively,  optical  fields  can  be  used  to  drive  two-photon  stimulated  Raman  transitions  between 
| n)  ||)j  and  | n')  ||)-  [17,  26].  As  depicted  in  Fig.  lb,  two  laser  beams  detuned  by  A  from  an 
excited  state  of  radiative  width  y  are  applied  to  the  jth  ion  with  their  difference  frequency 
matched  to  the  desired  transition  frequency.  For  sufficient  detuning  \A\  y,  the  excited 
state  may  be  adiabatically  eliminated,  and  the  above  couplings  apply,  with  Q!  replaced  by 
^1^2/ A,  where  gj  and  g,  arc  the  individual  Rabi  frequencies  of  the  two  beams  when  reso¬ 
nantly  coupled  to  the  excited  level.  In  addition,  m^cp)  is  replaced  by  the  difference  fre¬ 
quency  (phase)  of  the  beams;  and  k  is  replaced  by  the  difference  in  wavevectors  of  the  two 
Raman  beams  dk  =  k\  —  ki.  Since  the  relevant  frequency  depends  only  on  the  microwave 
difference  between  the  two  laser  frequencies,  both  beams  can  be  generated  with  a  single 
laser  source  and  a  modulator,  thereby  relaxing  the  constraints  of  laser  frequency  stabiliza¬ 
tion.  The  use  of  stimulated  Raman  transitions  thus  combines  the  advantages  of  the  strong 
couplings  with  the  frequency  stability  of  microwave  sources  [26,  27]. 


D.  Laser  cooling  to  the  motional  ground  state 

As  a  starting  point  for  trapped  ion  quantum  computing,  the  ions  must  be  initialized  in 
known  pure  states.  Using  standard  optical  pumping  techniques,  we  can  prepare  the  ions  in 
the  ||),  internal  state.  Laser  cooling  in  the  resolved  sideband  limit  [25]  can,  for  single  ions, 
generate  the  \n  =  0)  motional  state  with  reasonable  efficiency  [17,  28].  This  type  of  laser 
cooling  is  usually  preceded  by  a  stage  of  “Doppler”  laser  cooling  [29]  which  typically 
cools  the  ion  to  («)  >  1 ,  or  an  equivalent  temperature  of  about  1  mK. 

Resolved  sideband  laser  cooling  for  a  single,  harmonically-bound  atom  can  be  explained 
as  follows:  For  simplicity,  we  assume  the  atom  is  confined  by  a  1-D  harmonic  well  of 
vibration  frequency  atz.  We  use  an  optical  transition  whose  radiative  linewidth  y  is  rela¬ 
tively  narrow,  y  a>z  (Doppler  laser  cooling  applies  when  y  >  cu-).  If  a  laser  beam  (fre¬ 
quency  (>>[  )  is  incident  along  the  direction  of  the  atomic  motion,  the  bound  atom’s  absorp¬ 
tion  spectrum  is  composed  of  a  “carrier”  at  frequency  wq  and  resolved  frequency- 
modulation  sidebands  that  are  spaced  by  oaz.  These  sidebands  in  the  spectrum  are  generated 
from  the  Doppler  effect  (like  vibrational  substructure  in  a  molecular  optical  spectrum).  La¬ 
ser  cooling  can  occur  if  the  laser  is  tuned  to  a  lower  (red)  sideband,  for  example,  at 
a>i  =  o)q  —  a>z.  In  this  case,  photons  of  energy  h(a> 0  —  a>z )  are  absorbed,  and  spontaneously 
emitted  photons  of  average  energy  fia> 0  —  R  return  the  atom  to  its  initial  internal  state, 
where  R  =  (Tik)2 /2m  =  hwR  is  the  photon  recoil  energy  of  the  atom.  Overall,  for  each 
scattering  event,  this  reduces  the  atom’s  kinetic  energy  by  1ia>z  if  a>-  '»  coR,  a  condition 
which  is  satisfied  for  ions  in  strong  ion  traps.  Since  u>R/mz  =  rf  where  tj  is  the  Lamb- 
Dicke  parameter,  this  simple  form  of  sideband  cooling  requires  that  the  Lamb-Dicke  param¬ 
eter  be  small.  For  example,  in  9Be+,  if  the  recoil  corresponds  to  spontaneous  emission  from 
the  313  nm  2Pi/2  — >  2S!/2  transition  (typically  used  for  laser  cooling),  cor/2jt  ~  230  kHz. 
This  is  to  be  compared  to  trap  oscillation  frequencies  in  some  laser-cooling  experiments  of 
around  10  MHz  [17],  Cooling  proceeds  until  the  atom's  mean  vibrational  quantum  number 
in  the  harmonic  well  is  given  by  («)min  —  (y /2mff  1  [29].  As  discussed  above,  it  is 
convenient  to  use  two-photon  stimulated  Raman  transitions  for  sideband  cooling  [17,  30],  but 
the  basic  idea  for,  and  limits  to,  cooling  are  essentially  the  same  as  for  single-photon  transitions. 
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Although  laser  cooling  to  the  | n  =  0)  state  has  been  achieved  with  single  ions  [17,  28];  a 
prerequisite  to  future  work  is  to  laser  cool  a  collection  of  ions  (or,  at  least  one  mode  of  the 
collection)  to  the  zero-point  state.  Cooling  of  any  of  the  3 N  modes  of  motion  of  a  collec¬ 
tion  of  ions  should,  in  principle,  work  the  same  as  cooling  of  a  single  ion.  To  cool  a 
particular  mode,  we  tune  the  cooling  radiation  to  its  first  lower  sideband.  If  we  want  to 
cool  all  modes,  sideband  cooling  must  be  cycled  through  all  3 N  modes  more  than  once,  or 
applied  to  all  3 N  modes  at  once,  since  recoil  will  heat  all  modes.  For  the  COM  mode,  the 
cooling  is  essentially  the  same  as  cooling  a  single  particle  of  mass  Nm;  however,  the  recoil 
energy  upon  re-emission  is  distributed  over  all  3 N  modes.  Other  methods  to  prepare  atoms 
in  the  \n  =  0)  state  are  discussed  in  Refs.  [31].  In  Ref.  [32],  it  is  shown  that  it  is  not 
necessary  to  satisfy  the  condition  cdr  a>z  (rj  1)  to  achieve  cooling  to  n  =  0. 


III.  Quantum  Logic  with  Trapped  Ions 

Several  authors  have  shown  that  an  arbitrary  unitary  operation  (therefore  any  quantum  com¬ 
putation)  on  a  collection  of  quantum  bits  can  be  broken  into  a  series  of  fundamental  single 
bit  and  dual  bit  quantum  logic  gates  [33].  This  is  similar  to  classical  computing,  where 
certain  families  of  logic  gates  are  universal  (for  instance,  the  single  bit  NOT  and  two  bit 
AND  gates).  One  such  family  of  universal  quantum  logic  gates  consists  of  the  single  bit 
rotation  gate  and  the  two-bit  controlled-NOT  gate  [34].  For  brevity,  we  concentrate  on  these 
two  gates  and  how  they  can  be  implemented  in  a  system  of  trapped  ions. 

The  single  bit  rotation  gate  [operator  K(0.  cp)\  simply  changes  the  state  of  a  single  quan¬ 
tum  bit  and  is  characterized  by  the  following  transformation; 

11)^  cos  (0/2)  sin  (0/2)  It) 

IT)  ->  cos  (0/2)  It)  -  i  e~,rp  sin  (0/2)  | J.) , 

where  0  and  cp  are  parameters  of  the  gate.  This  transformation  is  commonplace  in  atomic 
physics  and  nuclear  magnetic  resonance,  and  has  been  widely  applied  to  two  level  systems. 
In  the  context  of  trapped  ions,  the  single  bit  rotation  gate  is  accomplished  by  tuning  to  the 
carrier  transition  ( a>i  =  coq)  and  applying  radiation  for  a  time  t  such  that  0  =  2f 2T 
(Eq.  (6)).  The  parameter  0  describes  the  “rotation"’  between  the  two  spin  states  (0  =  t  is 
called  a  At”  pulse,  etc.),  and  the  parameter  (p  describes  the  phase  of  the  rotation.  Figure  3 
depicts  the  observed  Rabi  flopping  between  ||),  and  |f),  states  in  a  single  trapped  9Be+  ion 


t(ns) 

Fig.  3:  Experimental  plot  of  the  probability  P\{t)  of  finding  a  single  9Be+  ion  in  the  ||) 
state  after  first  preparing  it  in  the  |J.)  state  and  applying  the  carrier  coupling  (Eq.  (7))  for  a 
time  t,  with  0/2  =  Q’t.  P[(t)  does  not  follow  a  perfect  sinusoid  due  to  decoherence,  de¬ 
scribed  in  section  IV.  Each  point  represents  an  average  of  4000  observations. 
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according  to  the  above  transformation  [35].  The  ion  is  first  prepared  in  the  |j),  state,  then 
the  carrier  transition  is  applied  for  a  time  t,  and  the  state  (||),  or  ||))  is  detected  as  de¬ 
scribed  above  and  in  Fig.  1  a. 

A  more  nontrivial  quantum  logic  gate  is  the  two-bit  controlled-NOT  (CN)  gate  [33,  34], 
reminiscent  of  the  classical  XOR  gate.  This  gate  flips  the  spin  of  ion  k,  the  “target  bit” 
(ID*  <->  |D*)  if  and  only  if  ion  j,  the  “control  bit,”  is  in  state  |]").  The  transformation  of  the 
two  quantum  bits  j  and  k  are  as  follows: 

\i)j  II}*  —  \i)j  II}* 

U-)j  IT)*  *•  IX)7-  IT)* 

IT);  ID*  —  It)*  IT)* 

IT}*  IT}*  -*•  |T)y  |i)jfc  - 

Cirac  and  Zoller  [9]  showed  how  the  above  transformation  could  be  accomplished 
between  two  of  a  collection  of  trapped  ions  by  utilizing  a  motional  mode  (i.e.,  axial  COM 
mode)  as  a  “data  bus”  through  which  the  quantum  bits  are  shuttled.  As  described  above, 
the  COM  mode  is  assumed  to  be  initially  cooled  to  the  | n  =  0)  ground  state.  The  scheme  is 
outlined  as  follows: 

(1)  Map  the  state  of  ion  j  (spanning  the  states  |J.) |]")p  onto  the  first  two  axial  COM 
motional  states  (spanning  the  states  |«  =  0),  |n  =  1))  shared  amongst  all  ions. 

(2)  Flip  the  spin  state  of  ion  k  (|D*  |T}*)  if  and  only  if  there  is  a  phonon  in  the  COM 

mode. 

(3)  Reverse  step  (1):  map  the  state  of  the  motion  back  onto  ion  j. 

The  central  ingredient  here  is  step  (2),  which  is  itself  a  “reduced”  CN  gate,  with  the  first 
two  motional  states  acting  as  the  control  bit  and  ion  k  acting  as  the  target  bit 

|o)  ID* -10}  ID* 

10}  IT}*  -  10}  IT}* 

|i}  II}* -|i}  IT}*  1  j 

|i}  IT}*—  |i}  ID*- 

This  last  transformation  has  been  realized  on  a  single  trapped  ion  [10].  In  that  experiment, 
performed  on  a  trapped  9Be+  ion,  the  control  bit  was  one  of  the  three  modes  of  the  ion’s 
motion.  The  reduced  CN  operation  between  these  states  (step  (2)  above)  was  realized  by 
applying  three  laser  pulses  in  succession: 

(a)  A  n j 2  pulse  ( Qk  t  =  .t/4  in  Eq.  (6))  is  applied  on  the  carrier  transition.  For  a  certain 
choice  of  initial  phase,  this  corresponds  to  the  rotation  operator  R(6  =  ti/2,  <p). 

(b)  A  2jt  pulse  is  applied  on  the  first  blue  sideband  transition  between  levels  ||}  and  an 
auxiliary  level  |aux)  in  the  ion  (the  \F  =  2,  Mp  =  0}  level  in  9Be+;  see  Fig.  4).  This 
operation  provides  the  “conditional  dynamics”  for  the  overall  CN  operation.  It  changes 
the  sign  of  the  |1)  |j)  component  of  the  wavefunction  but  leaves  the  sign  of  the  |0)  |T} 
component  of  the  wavefunction  unchanged;  that  is,  the  sign  change  is  conditioned  on 
whether  or  not  the  ion  is  in  the  |0)  or  |l)  motional  state. 

(c)  A  jr/2  pulse  is  applied  to  the  spin  carrier  transition  with  a  180°  phase  shift  relative  to 
step  (a).  This  corresponds  to  the  operator  R(8  =  ti/2,  cp  +  ji). 

Steps  (a)  and  (c)  can  be  regarded  as  two  resonant  pulses  of  opposite  phase  in  the  Ramsey 
separated-field  method  of  spectroscopy  [36].  If  step  (b)  is  active  (thereby  changing  the  sign 
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(a)  (b)  (c) 

Fig.  4:  Energy  levels  of  a  single  trapped  9Be+  ion,  including  the  internal  hyperfine  levels 
||)  and  ||)  and  an  auxiliary  level  |aux)  [the  2Si/2(F  =  2,  iiif  =  0)  state],  each  dressed  by 
the  lowest  two  motional  quantum  harmonic  oscillator  states  | n)  =  |0)  and  |1).  The  con- 
trolled-NOT  quantum  logic  gate  results  in  a  spin  flip  (||)  ||)  if  and  only  if  | n)  =  |1). 

This  transformation  is  realized  with  a  sequence  of  three  pulses  of  laser  light  which  couple 
the  states  indicated  by  the  arrows:  (a)  A  jr/2  pulse  couples  states  |/i)  ||)  to  | n)  ||).  (b)  A  2 n 
pulse  couples  state  1 )  ||)  to  |0)  |aux),  resulting  in  a  sign  change  of  any  component  in  the 
|1)  |t)  state,  (c)  A  —n/2  pulse  couples  states  \n)  ||)  to  \n)  ||)  (same  as  step  (a)  with  a  jt 
phase  shift).  If  \n)  =  |0),  then  step  (b)  is  inactive  since  it  only  affects  to  the  |1)  ||)  state,  and 
the  two  jr/2  pulses  cancel,  leaving  the  initial  state  unaffected.  If  |/t)  =  |l),  then  the  sign 
change  in  step  (b)  causes  the  two  jr/2  pulses  to  add,  resulting  in  a  net  spin  flip  (||)  ±=>  ||)). 
The  result  is  the  transformation  of  Eq.  (9). 


of  the  |1)  ||)  component  of  the  wavefunction),  then  a  spin  flip  is  induced  by  the  Ramsey 
fields.  If  step  (b)  is  inactive,  step  (c)  reverses  the  effect  of  step  (a). 

Instead  of  the  three  pulses  (a— c  above),  a  simpler  CN  gate  scheme  between  an  ion’s 
internal  and  motional  states  can  be  achieved  with  a  single  laser  pulse,  while  eliminating  the 
requirement  of  the  auxiliary  internal  electronic  level  [37].  By  applying  a  single  pulse  tuned 
to  the  carrier  transition,  from  Eq.  (5),  the  states  | n)  ||)  and  | n)  ||)  are  coupled  with  Rabi 
frequency 

Q„,n  =  &  IH  ^(fl+flt)  |n)|  =  Qe-*2'2  J?n(V2) ,  (10) 


where  £> ?n(rj2)  =  L®(rj2)  is  a  Laguerre  polynomial.  Specializing  to  the  \n)  =  |0)  and 
| n)  =  |1)  vibrational  levels  relevant  to  quantum  logic,  we  have 


Q0,o  =  tie-'12 12' 
f 2,.!  =  Q e~'l2!2{\  —  rj2) . 


(11) 


The  CN  gate  can  be  achieved  in  a  single  pulse  by  setting  77  so  that 
®i,i/^o,o  =  {2k  +  l)/2m,  with  k  and  m  positive  integers  satisfying  m>k>  0.  Setting 
£2i,i/£2o.o  =  2m/(2k  +  1)  will  also  work,  with  the  roles  of  the  |0)  and  |1)  motional  states 
switched  in  Eq.  (10).  By  driving  the  carrier  transition  for  a  duration  r  such  that 
iT  =  (k  + Jt,  or  a  “jr-pulse”  (mod  2jt)  on  the  \n)  =  |1)  component,  this  forces 
£2q  qZ  =  rrm.  Thus  the  states  ||)  |1)  |T)  |1)  are  swapped,  while  the  states  ||)  |0)  and 

|t)  |0)  remain  unaffected.  This  transformation  is  equivalent  to  the  reduced  CN  of  Eq.  (9), 
apart  from  phase  factors  which  can  be  eliminated  by  the  appropriate  settings  of  the  phase 
of  subsequent  logic  operations. 

The  mapping  steps  (steps  (1)  and  (3)  above)  can  be  realized  by  applying  a  jr-pulse  on 
the  red  sideband  of  ion  j.  This  accomplishes  the  mapping  |0)  (a  ||)  .  +  /?  |t)y-)  — >  (a  |0) 
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+  /3|1))  Analogous  mapping  of  internal  state  superpositions  to  motional  state  super¬ 
positions  of  a  trapped  ion  has  also  been  reported  in  the  generation  of  a  “Schrodinger 
cat”  state  of  motion  [38]  and  the  tomographic  measurement  of  nonclassical  states  of 
motion  [39]. 

To  complete  this  section,  we  give  an  example  of  how  a  “maximally-entangled"  state  of  N 
ions  might  be  prepared  with  the  use  of  the  quantum  logic  gates  described  in  this  section 
[9].  We  desire  to  create  the  state 


^*ME  — 


li)i  li)2  U>3  -  -  -  U)at  +  ^  IT)t  lt)2  lt)3  -  -  •  IT)at 
V2 


(12) 


(or  equivalently,  a  coherent  superposition  of  the  numbers  0  and  2N  —  1  in  binary  where  we 
make  the  identification  ||)  =  0  and  |j)  =  1).  Starting  in  an  initial  state  0init  where  all  spins 
are  in  the  |j)  state,  it  is  easy  to  show  that  one  way  to  prepare  Eq.  (12)  is  to  operate  the 
following  N  gates  on  the  initial  state : 


<2>me=  [CM,iv][CAfiijv-1][CM,Ar-2]...[CiV1,2][/?i(jr/2)  q>)\0 m,  (13) 


where  CNjj  denotes  the  controlled-NOT  operator  with  ion  i  as  the  control  bit  and  ion  j  as 
the  target  bit,  and  the  rotation  operator  is  applied  to  the  first  ion.  The  state  <Pme  can  be 
viewed  as  the  /V-particle  generalization  of  the  entangled  pair  states  envisioned  by  Einstein, 
Podolsky,  and  Rosen  [40],  and  may  find  uses  in  improved  atomic  spectroscopy  and  fre¬ 
quency  standards  [7,  8]. 


IV.  Packing  Ions  into  a  Trap 

The  reduced  CN  gate  demonstrated  in  [10]  involved  only  a  single  ion,  and  is  therefore  not 
useful  for  computation.  Future  experimental  work  will  concentrate  on  scaling  the  system  up 
by  packing  more  ions  into  the  trap  and  improving  the  gate  fidelity,  thereby  allowing  more 
gates  to  be  coherently  performed.  In  this  section,  we  survey  selected  experimental  problems 
which  may  arise  in  the  scaleup.  For  a  more  extensive  coverage  of  the  scaling  problem,  see 
Ref.  [11], 


A.  Individual  ion  addressing 


One  major  concern  in  extending  quantum  logic  to  larger  numbers  of  ions  is  the  requirement 
that  ions  be  individually  addressed  with  laser  beams  for  logic  operations.  This  may  be 
difficult,  because  the  high  vibrational  frequencies  desired  for  efficient  laser  cooling  and 
suppression  of  decoherence  also  results  in  closely  spaced  ions.  It  can  be  shown  that  the 
minimum  separation  of  adjacent  ions  in  a  linear  trap  between  the  center  ions  scales  approxi¬ 
mately  as  vmin  ~  2sIV_0'  ,  with  s  =  (q2 /4ji£oma>l)1'3 ,  where  N  is  the  number  of  ions 
[12,  13].  For  9Be+  ions  with  an  axial  COM  frequency  of  cjz/2jz  =  1  MHz,  this  separation 
is  about  10  pm  for  2  ions,  and  4  pm  for  10  ions. 

The  most  straightforward  method  for  individual  optical  addressing  is  to  tightly  focus 
laser  beams  on  the  selected  ion  [9].  The  transverse  intensity  distribution  of  a  Gaussian 
optical  beam  of  power  P  is 


Hr)  = 


2  P 
nw. 


exp 


2  r2 

wl 


(14) 
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where  wq  —  •  NA)  is  the  beam  waist,  A  the  radiation  wavelength,  and  NA  =  tan  0  is  the 

numerical  aperture  of  the  beam  with  cone  half-angle  6  (the  formula  for  h’o  in  the  paraxial 
ray  approximation  is  valid  only  for  NA  <  1)  [41].  For  large  numerical  apertures 
(NA  «  0.5),  it  appears  that  laser  beams  can  thus  be  focussed  down  to  a  spot  on  the  order  of 
a  wavelength,  but  this  is  difficult  to  realize  in  the  laboratory.  If  we  can  realize  wo  =  5  pm 
in  a  Gaussian  beam,  then  at  a  distance  10  pm  from  the  center  of  the  beam  this  would 
imply  a  relative  intensity  of  about  3  •  1 0  4  or  a  electric  field  amplitude  (proportional  to 
Rabi  frequency)  of  1.8%  relative  to  the  center  of  the  beam.  If  wo  =  2  pm  could  be  ob¬ 
tained,  the  intensity  (electric  field)  would  be  down  by  a  factor  of  1.3  •  10~14  (1.1  •  10  7). 
These  results  are  likely  too  optimistic,  since  imperfections  in  the  surfaces  of  the  intervening 
vacuum  port  window,  multiple  reflections  from  these  windows,  and  diffraction  will  typi¬ 
cally  distribute  laser  intensity  outside  of  the  theoretical  waist  of  the  beam.  The  degree  to 
which  this  occurs  depends  on  the  details  of  window  surfaces,  etc.  and  must  be  resolved 
experimentally. 

Tightly  focussed  Gaussian  beams  possess  high  transverse  intensity  gradients,  resulting  in 
the  potential  for  significant  intensity  fluctuations  at  the  selected  ion  if  the  relative  position 
of  the  beam  with  respect  to  the  ion  is  not  stable  on  the  time  scale  of  the  computation.  An 
alternative  to  using  tightly  focussed  Gaussian  laser  beams  is  to  first  feed  the  (expanded) 
laser  beam  through  a  sharply  defined  aperture  (slit  or  aperture),  and  use  a  lens  to  image  the 
aperture  at  the  position  of  the  ions.  With  this  technique,  the  beam  intensity  can  be  distribu¬ 
ted  more  smoothly  around  the  selected  ion  and  have  very  steep  intensity  edges  (on  the 
order  of  the  original  aperture  sharpness)  away  from  the  ion,  thus  suppressing  beam  vibra¬ 
tion  problems  and  confining  the  radiation  to  a  single  ion.  This  technique  has  been  used  to 
make  relatively  “hard”  walls  for  an  optical  dipole  trap  [42].  For  this  technique  to  work 
well,  the  imaging  lens  must  collect  a  large  fraction  of  the  light  transmitted  through  the 
aperture  or  else  diffraction  effects  will  result  in  light  intensity  outside  the  image  of  the 
aperture.  To  address  individual  ions,  we  require  very  small  aperture  images,  which  gives 
rise  to  a  design  tradeoff.  If  a  one-to-one  relay  lens  is  used  to  image  a  small  object  aperture, 
effects  of  diffraction  are  enhanced.  If  a  demagnifying  lens  is  used  to  reduce  a  large  object 
aperture,  then  the  aperture  must  be  placed  a  large  distance  from  the  lens,  requiring  a  rela¬ 
tively  large  lens.  For  two  ions,  imaging  a  sharp  edge  such  as  a  razor  blade  at  the  space 
between  the  ions  may  be  sufficient.  We  might  also  consider  having  every  other  ion  in  a 
string  be  a  “garbage”  ion  which  is  not  used  in  the  computation,  thereby  increasing  the 
spacing  between  qubit  ions  by  a  factor  of  two  (or  more,  if  more  garbage  ions  are  used 
between  each  qubit  ion).  This  has  the  disadvantage  that  total  number  of  ions  (and  spectator 
modes)  increases,  aggravating  the  problems  associated  with  large  quantum  registers.  If  suffi¬ 
ciently  good  addressing  on  one  ion  in  a  string  can  be  accomplished,  it  may  be  simpler  to 
shift  the  ions,  rather  than  the  laser  beams,  in  order  to  address  different  ions.  This  could  be 
accomplished  by  applying  different  static  potentials  Uq  and  U'Q  to  the  end  segments  of  the 
rods  in  Fig.  2.  However,  changes  in  Uq  and  Ur0  would  have  to  be  coordinated  to  keep  the 
COM  axial  frequency  constant  or  else  additional  phase  shifts  would  be  introduced.  Stimu- 
lated-Raman  transitions  have  the  advantage  that  the  effective  wavevector  dk  =  k\  —  k2  can 
be  made  parallel  to  the  axis  of  the  trap  even  though  each  beam  is  at  an  angle  with  respect 
to  the  trap  axis.  This  would  allow  selection  of  a  particular  ion,  while  eliminating  coupling 
to  transverse  modes. 

Another  method  of  optically  addressing  individual  ions  is  to  cause  a  destructive  optical 
interference  at  the  position  of  a  specific  ion,  with  a  net  coupling  at  the  other  ion(s).  For 
instance,  if  ion  j  is  positioned  at  the  node  of  a  resonant  standing  wave  laser  field,  the 
coupling  between  states  | n)  ||)  -  and  | n')  |"f)-  is  proportional  to  (n'\  sin  [r]j(a  +  a')]  | n).  In  this 
case,  the  coupling  of  the  standing  wave  to  ion  j  vanishes  when  the  laser  frequency  is  tuned 
to  an  even  order  sideband  such  as  the  carrier  («'  =  n).  If,  instead,  the  ion  is  positioned  at 
an  antinode,  the  coupling  is  proportional  to  (n'|  cos  [r]j(a  +  a')]  |n);  thus,  the  coupling 
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vanishes  when  the  laser  frequency  is  tuned  to  an  odd  order  sideband,  such  as  the  first  blue 
or  red  sideband  («'  =  n±  1).  By  choosing  the  angles  of  focussed  laser  beams  relative  to 
the  trap  axis  appropriately,  it  should  be  possible  to  position  an  antinode  (node)  at  ion  j 
while  approximately  positioning  nodes  (antinodes)  at  the  ions  adjacent  to  ion  j  (for  equally 
spaced  ions).  In  the  case  of  two-photon  stimulated  Raman  transitions,  we  desire  to  place 
ion  j  at  a  common  node  or  antinode  of  two  standing  waves.  Although  this  interference 
technique  should  allow  individual  access  to  each  of  two  trapped  ions,  it  appears  technically 
difficult  to  extend  this  technique  to  more  than  three  ions. 

Next,  we  consider  the  application  of  external  field  gradients  which  shift  the  internal 
energy  levels  of  ions  depending  on  their  position.  For  a  magnetic  field  gradient  to  give  this 
selectivity,  we  require  the  Zeeman  splitting  between  adjacent  ions  to  be  much  larger  than 
the  Rabi  frequency,  or  Afi{d\B\/dz)  s/h^>  Q,  where  A/i  is  the  difference  in  {[i-B)/\B\ 
between  the  two  levels  of  interest,  and  s  is  the  ion  separation  along  the  z  direction.  For 
Afi  =  jUB ,  s  =  1 0  um.  and  a  Rabi  frequency  of  Q lire  I  MHz,  this  requires  <9|/?|/<9z 

0.1  T/cm.  Field  gradients  of  this  magnitude  can  be  achieved;  however,  they  might  intro¬ 
duce  large  and  uncontrollable  phase  shifts  for  the  other  ions  in  a  quantum  register. 

The  laser  beam  itself  can  provide  ion  selectivity  by  employing  the  transverse  gradient  in 
the  optical  field  intensity.  For  instance,  if  we  desire  to  perform  a  0-pulse  on  ion  j  without 
affecting  neighboring  ion  k,  the  intensity  profile  of  the  laser  beam  can  be  set  so  that  the 
ratio  of  field  strengths  (intensities  for  the  case  of  two-photon  stimulated-Raman  transitions) 
at  ion  j  vs.  ion  k  is  6 /2mn,  where  m  is  an  integer.  Now  if  the  pulse  duration  is  set  so  that 
ion  j  is  rotated  by  0,  ion  k  receives  a  rotation  of  2mn  and  hence  returns  back  to  its  initial 
state  (with  an  extra  phase  factor  of  (—1)"'). 

For  the  case  of  two-photon  stimulated-Raman  transitions,  the  laser  beam  can  provide  ion 
frequency  selectivity  by  employing  the  Stark  shift  and  the  transverse  gradient  of  the  optical 
field.  Here,  for  example,  we  could  assume  that  the  two  counterpropagating  Raman  beams 
of  equal  intensities  and  spatial  profiles  are  offset  so  that  beam  1  is  centered  on  ion  j ,  and 
beam  2  is  centered  on  adjacent  ion  k  as  depicted  in  Fig.  5.  Let  e  be  the  fraction  of  peak 
intensity  seen  by  the  offset  ions  (that  is,  the  intensity  of  beam  2  at  ion  j  and  beam  1  at  ion 
k).  Assume  that  when  either  beam  is  centered  on  either  ion,  the  single  photon  resonant 
Rabi  frequencies  are  equivalent:  g\  =  g2  =  g-  When  the  beams  are  offset,  the  two-photon 
Rabi  frequency  at  each  ion  is  Q  =  £l'2(g2  / Ar),  where  g2 /Ar  is  the  Rabi  frequency  ex¬ 
pected  if  both  beams  were  centered  on  a  given  ion.  The  Stark  shifts  of  the  two  ions  are  in 
opposite  directions:  6 j  =  +(5o,  <5^  =  — <5o-  where  do  =  Q{\  —  ej/e1/2.  If  we  make  Q 

(e  1),  then  by  appropriately  tuning  the  difference  frequency  of  the  laser  beams,  we 
can  selectively  drive  transitions  on  either  ion  j  or  k.  Alternatively,  if  we  desire  to  per¬ 
form  a  0-pulse  to  ion  j  without  affecting  ion  k  in  an  “unrepairable”  way,  £  can  be 
tuned  to  a  particular  value  which  results  in  ion  k  returning  to  within  a  phase  factor  to 
its  initial  state.  For  square  pulses  in  time,  we  require  0/2(1  +  d2/Q~)l^~  =  rrm,  or 
£2  —  [1  +  (2mJi/6)2]  £+1=0,  where  m  is  an  integer.  For  m  =  1  and  0  =  j r  (a  jr-pulse  on 
ion  j),  this  occurs  for  £  =  0.208.  Generalizing  this  to  more  than  two  ions  becomes  difficult 
if  the  laser  beams  also  overlap  other  qubit  ions.  This  scheme  places  an  additional  premium 
on  laser  power  stability,  since  the  light  shifts  are  bigger  than  the  Rabi  frequencies  by  1/e1/2 
for  £  1.  In  addition,  in  both  of  the  above  schemes,  employing  the  laser  beams  to  differ¬ 

entially  affect  neighboring  ions,  one  major  drawback  is  that  the  positions  and  profiles  of  the 
laser  beams  must  be  accurately  controlled. 

Finally,  we  consider  a  method  of  addressing  which  utilizes  rf  micromotion.  In  a  con¬ 
ventional  Paul  trap,  generalized  to  the  case  of  asymmetric  electrodes  (0  oc  Vo  cos  (Qjt) 
x(ex2  +  (1  —  £)y2  —  z2),  0  <  £  <  1)  [43],  the  rf  fields  vanish  only  at  a  single  point  in 
space.  When  multiple  ions  are  crystalized  in  such  a  trap,  each  ion  experiences  a  different  rf 
field,  leading  to  different  amounts  of  micromotion.  In  general,  this  effect  causes  each  ion  to 
have  a  unique  Rabi  frequency  when  a  laser  is  applied,  allowing  the  possibility  of  differen- 
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Fig.  5 :  Differential  AC  Stark  shift¬ 
ing  of  neighboring  ions.  Equal  in¬ 
tensity  counter  propagating  beams  1 
and  2  are  centered  on  ions  j  and  k, 
respectively.  A  fraction  e  of  the 
peak  intensity  /o  of  each  beam  is 
applied  to  the  other  ion.  This  results 
in  a  differential  AC  Stark  shift  of 
ions  j  and  k ,  allowing  the  possibility 
of  individually  accessing  the  ions  by 
tuning  the  frequency  of  the  laser 
beams. 


tial  addressing  of  the  ions.  If  we  assume  that  each  ion  sees  the  same  intensity  laser  field, 
the  Rabi  frequency  of  ion  j  is  proportional  to  Ja(k  ■  £)■)  [1 1,  44],  where  Jo(X)  is  the  zeroth 
Bessel  function  of  argument  X,  k  is  the  effective  laser  wavevector,  and  is  the  amplitude 
of  micromotion  associated  with  the  /th  ion.  We  assume  that  ion  j  experiences  a  static  elec¬ 
tric  field  Esj  (due  to  the  Coulomb  field  of  the  other  ions  in  addition  to  background  static 
fields)  which  prevents  the  ion  from  occupying  the  trap  center  defined  by  the  rf  fields.  For 
simplicity,  we  assume  that  Esj  =  ESjX  is  along  the  x  axis  of  the  trap,  which  should  be  the 
case  if  the  ions  crystalize  along  the  x-axis  (e  <  0.5).  To  estimate  u,  we  solve  the  classical 
equations  of  motion  in  the  trap  for  Uo  =  0  and  find  =  ^  •  x  =  \/2qESj/ma>xQT,  where  mx 
is  the  x-axis  COM  secular  harmonic  frequency.  Thus,  the  Rabi  frequency  of  the  /th  ion  is 


Qj  =  Q'()J0 


(15) 


where  Q'(]  is  the  Rabi  frequency  in  absence  of  micromotion.  (This  reduction  of  the  Rabi 
frequency  due  to  ion  motion  is  treated  in  section  IV.B.l.  under  a  different  context.)  Thus  by 
controlling  the  electric  field  Ej  at  the  /th  ion,  some  degree  of  differential  addressing  of  the 
ions  is  possible.  For  example,  for  three  ions  held  along  the  x-axis  with  the  middle  ion 
placed  exactly  at  the  rf  null  position,  we  find  that  the  ratio  of  Rabi  frequencies  of  the 
outer  ions  to  the  middle  ion  is  Jq(/3),  where  /3  =  l.52kx(q2cox/4jt£om)1^/ Qj.  An  interest¬ 
ing  case  for  individual  addressing  occurs  when  /3  =  2.405,  in  which  case  Jo(fi)  =  0,  and 
the  laser  interaction  with  the  outer  two  ions  is  effectively  shut  off.  For  three  9Be+  ions 
held  in  a  trap  with  a  rf  drive  frequency  of  Qj  «  240  MHz,  we  find  that  the  Rabi  frequency  of 
the  outer  two  ions  vanishes  at  a  COM  secular  frequency  of  cox/ 2?r  «  6.1  MHz.  Here  we 
assume  that  the  Rabi  frequency  describes  a  stimulated  Raman  coupling  between  hyperfine 
ground  states  with  kx  =  \/2(2n/X)  and  A  «  313  nm,  relevant  to  recent  experiments  [10]. 

Many  of  the  above  individual  addressing  schemes  are  improved  greatly  when  dealing 
with  only  two  ions  instead  of  a  string  of  many.  This  leads  us  to  seriously  consider  systems 
where  quantum  logic  operations  are  performed  on  accumulators  consisting  of  only  two 
ions,  with  the  other  ions  located  somewhere  else.  For  example,  pairs  of  ions  may  be  held  in 
different  regions  of  the  same  trap  structure  [11],  or  quantum  information  may  be  transferred 
from  one  register  of  ions  to  another  by  optical  means  [45]. 
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B.  Multimode  interference 


Each  additional  ion  in  a  quantum  register  adds  three  more  motional  modes.  The  simplest 
form  of  the  Cirac/Zoller  scheme  ideally  uses  just  a  single  mode  as  the  bus.  In  this  section, 
a  few  potential  problems  associated  with  the  3 N  —  1  spectator  modes  are  considered.  First, 
we  must  generalize  the  interaction  with  electromagnetic  fields  discussed  in  section  II. B  to 
consider  motion  in  all  3 N  modes  of  motion  for  N  trapped  ions.  Here,  as  was  assumed  by 
Cirac  and  Zoller  [9],  we  consider  that,  on  any  given  operation,  the  laser  beam(s)  interacts 
with  only  the  /'th  ion;  however,  that  ion  will,  in  general,  have  components  of  motion  from 
all  modes.  In  this  case  Eq.  (3)  for  the  /th  ion  becomes 

H\  =  TiQ\S\  +  Sj_)  [eK*-x-<ot-Hpj)  +  h  c  ]  _  (16) 


We  write  the  position  operator  of  the  /'th  ion  (which  represents  the  deviation  from  its  equili¬ 
brium  position)  as 

Xj  =  ujx +  uN+jy +  U2N+jZ,  j&  {1,  2, ,...  N},  (17) 

where  the  up  are  related  to  the  3N  normal  mode  coordinates  qk  (k  €  {1,  2, . . .  3N}  by  the 
following  relations  [46] 


Up  =  Y  DkVk  ,  qk=  Y  Hfap  >  <lk  =  qm{ak  e  ,0Jk‘  +  a\  e,a,t') .  (18) 

k—\  p= 1 

In  this  expression,  qk  is  the  position  operator  and  ak  and  a\  are  the  associated  lowering  and 
raising  operators  for  for  the  A:th  normal  mode,  and  the  matrix  ifk  is  the  transformation 
matrix  between  physical  coordinates  of  the  individual  ions  and  normal  coordinates  of  the 
string.  Following  the  procedure  of  section  II. C,  we  find 


H\  =  hQ'  exp 


i  Y  rft (ake  KOk‘  +  al  e,a>kt)  —  i(dt  —  cpj)  +  h.c.  ]  , 

.  *  =  l  J  J 


(19) 


where  we  have  kept  the  time  dependence  to  allow  consideration  of  all  motional  modes  with 
different  frequencies,  and  d  =  coL  —  a>0  is  the  detuning  of  the  applied  radiation  frequency 
(or  difference  frequency  for  the  Raman  coupling).  The  generalized  Lamb-Dicke  parameters 
are  rfk  =  (k  ■  xiyk  +  k  ■  yD^+j  +  k  ■  zD2kN+j)  quo,  but  for  the  linear  trap  case,  motion  will  be 
separable  in  the  x,  y,  and  z  directions  and  rfk  will  consist  of  one  of  these  terms.  We  are 
typically  interested  in  coupling  the  internal  states  of  a  given  ion  j  to  a  selected  mode  of  collec¬ 
tive  motion  k.  In  this  case,  the  Rabi  frequency  coupling  the  states  \nk)  |J.)  -  to  | n'k)  |t);-  is 


,  3  N  ,  t 

nk  =  & nk I  n  ^'(a,+a']  I {nptkh  nk) I  • 
k’  1=1 


(20) 


where  \{np^k})  denotes  the  state  of  motion  of  the  2>N  —  1  spectator  modes  excluding 
mode  k. 


1.  Effects  of  motion  in  spectator  modes  on  logic  gates  (Debye- Waller  factors) 

From  the  last  equation,  the  Rabi  frequency  of  a  particular  operation  will  in  general  depend 
on  the  motional  state  of  the  spectator  modes.  For  instance,  the  conventional  controlled-not 
gate  employs  two  carrier  pulses  (steps  (a)  and  (c)  in  Sec.  Ill)  which  ideally  do  not  depend 
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on  the  state  of  motion;  this  requires  the  Lamb-Dicke  parameter  i)  to  be  small  (see  Eq.  (10)). 
In  the  Raman  configuration,  r\  is  proportional  to  the  difference  in  two  wavevectors  and  can 
be  made  negligible  by  using  co-propagating  beams  (dk  w  0).  On  the  other  hand,  with  sin¬ 
gle-photon  optical  transitions,  the  Rabi  frequencies  depend  on  the  motion  of  all  modes 
which  have  a  component  of  motion  along  the  direction  of  k.  We  can  take  advantage  of  the 
motional  dependence  of  the  carrier  to  construct  a  logic  gate,  but  in  this  case  also,  the  Rabi 
frequency  will  depend  on  the  motion  in  the  other  modes  along  the  direction  of  k  or  dk. 
Similarly,  for  sideband  operations,  such  as  steps  (1)  and  (3)  of  the  CN  scheme  discussed  in 
Sec.  Ill,  it  will,  in  general,  be  impossible  to  have  Q'n,  depend  on  only  one  mode  of 
motion.  In  this  section,  we  examine  the  influence  of  extraneous  modes  on  the  Rabi  frequen¬ 
cies  Q'n,n. 

In  a  collection  of  N  ions,  the  motion  in  the  3 N  —  1  spectator  modes  reduces  the  Rabi 
frequency  in  much  the  same  way  as  lattice  vibrations  affect  a  single  emitter  or  scatterer 
embedded  in  a  crystal,  as  described  by  the  Debye-Waller  effect  [47].  Typically,  the  mo¬ 
tional  quantum  numbers  of  the  spectator  modes  in  Eq.  (20)  will  be  thermally  distributed 
(i.e.,  P(iip)  oc  yn,‘  where  y  =  np/(  1  +  hp)  with  np  the  mean  number  of  phonons  in  mode  p), 
so  we  can  calculate  the  rms  and  mean  values  of  the  Rabi  frequency  given  this  distribution. 
For  simplicity,  we  assume  that  all  spectator  modes  are  in  the  Lamb-Dicke  regime 
{rf  nlJ2  «  1),  but  see  Ref.  [11]  for  the  more  general  case.  If  the  frequencies  and  ampli¬ 
tudes  of  all  modes  contributing  to  the  axial  motion  of  ion  j  are  assumed  to  be  about  the 
same,  we  can  write  rfp  ~  ij/N1'2  and  np  ~  n  where  ij  and  h  are  the  Lamb-Dicke  parameter 
and  mean  occupation  for  the  axial  motion  of  a  single  (thermalized)  trapped  ion.  In  this 
case,  the  fractional  fluctuation  in  the  Rabi  frequency  from  run  to  run  is 


aq1^ 

n>k  ~  rP- 

O’, 


n(n  +  1) 


(21) 


This  expression  indicates  that  a  large  number  of  ions  is  beneficial  because  it  tends  to  aver¬ 
age  out  the  effects  of  motion  in  the  N  —  1  spectator  modes.  Eq.  (21)  is  an  overestimate  of 
the  fluctuations  since  the  spectator  modes  will  have  higher  frequency  than  the  COM  mode, 
leading  to  smaller  amplitudes  of  motion  than  assumed  in  this  crude  estimate.  In  any  case,  it 
is  clearly  desirable  to  cool  all  modes  (whose  motion  is  parallel  to  dk)  to  the  zero-point 
state  0ip  =  0)  to  suppress  the  effects  of  these  Debye-Waller  factor  fluctuations. 


2.  Mode  cross-coupling  from  static  electric  field  imperfections 

If  the  3 N  —  1  spectator  modes  of  oscillation  are  not  all  laser-cooled  to  their  zero-point 
energy,  then  energy  can  be  transferred  to  the  kth  mode  of  interest.  Even  when  the  spectator 
modes  are  cooled  to  the  zero-point  state,  they  can  act  as  a  reservoir  for  energy  from  the 
COM  mode.  This  does  not  lead  to  heating  but  can  cause  decoherence. 

Ideally,  the  ions  are  subjected  to  quadratic  potentials  as  in  Sec.  II.B.  In  practice,  higher- 
order  static  potential  terms  are  present;  these  terms  can  induce  a  coupling  between  the 
modes.  Similar  couplings  are  induced  by  the  time  varying  fields  necessary  for  providing 
entanglement;  these  are  discussed  below.  We  will  assume  that  the  higher  order  field  gradi¬ 
ents  act  as  a  perturbation  to  the  (harmonic)  normal  mode  solution.  Following  the  conven¬ 
tion  in  the  above  introduction  to  Sec.  IV.B,  these  fields  will  be  specified  by  E,  for 
i  €  {1,  2,  ...  3N}  where  the  index  i  specifies  both  the  ion  and  direction  of  E.  We  write  the 
electric  field  at  the  /th  ion  as 


Ej  —  EjX  +  EN+jy  +  E2N+jz ,  j  €  {1,  2,  . . .  N}  . 


(22) 
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From  Eqs.  (18)  and  (21),  we  can  write  the  equation  of  the  Ath  normal  mode  as  [46] 

3  N 

E 

f=l 


(23) 


In  general,  we  can  write 


Ei  =  mup })  =  Ei({qj}) 

3  N 

=  Ei{{qj}  =  0)  +  J2 

m  —  1 


dEi 


dqIT 


-I  UA  =  o 


l  3  L  3N 

-W  E  E  mrr 
^  /—  1  m  —  1 


92£; 


%  %> 


J  {^}  =  0 


(24) 


where  the  derivatives  are  evaluated  at  the  equilibrium  positions.  The  first  term  on  the  right 
side  of  the  last  equation  just  gives  rise  to  a  shift  of  the  equilibrium  positions,  and  the 
second  term  can  be  absorbed  into  new  definitions  of  the  normal  mode  frequencies  «),.  The 
second-order  term  (last  term  shown  in  this  equation)  can  resonantly  couple  two  modes  of 
oscillation  (Z  and  m)  to  the  normal  mode  of  interest  k.  We  find  a  possible  resonant  term: 


d2qk 

df- 


+  <■ t>~kqk 


q  3N 

E  D/.q/qm 

YYl  i  —  i 


d2E, 


dqi  dqm\  {?.}  =  0 


(25) 


where  the  Z  and  m  mode  frequencies  satisfy  \a>i  ±  com\  =  a>k.  This  type  of  coupling 
can  either  add  to  or  extract  energy  from  mode  k,  depending  on  the  relative  phases 
of  motion  in  the  three  modes.  By  substituting  the  free  solution  to  modes  Z  and 
m  (qj(t)  =  Qk  exp  (±i(a)jt  +  qpj))  into  the  second-order  term,  we  find  that  if 
qk  (t  =  0)  =  (dqk/dt)t=0  =  0,  the  driven  solution  to  the  amplitude  of  mode  k  initially 
grows  linearly  with  time: 


qt 

2  ma>k 


3  N 

E  D'kQ,Qm 


i=  1 


d2Ei 

dq,  dqm 

{«}  =  ° 

(26) 


We  illustrate  with  an  approximate  numerical  example  which  might  have  been  expected  to 
play  a  role  in  the  motional  heating  that  was  observed  in  the  NIST  experiments  [10,  17] 
(see  section  V.B  below).  In  those  experiments,  performed  on  single  4Be+  ions,  the  heat¬ 
ing  that  was  observed  on  the  motion  in  the  x  direction  was  such  that  the  ion  made  a 
transition  from  the  n  =  0  to  n  =  1  level  in  about  1  ms.  For  a  single  ion,  the  three  nor¬ 
mal  modes  are  just  the  oscillation  modes  along  the  x,  y,  and  z  directions  (q i  =  x, 
q2=y,  qi  =  z;  D'k  =  dk).  The  mode  frequencies  were  (cyx,  a>y,  (dz)/2jt  ~(11.2,  18.2, 
29.8)  MHz,  thus  approximately  satisfying  the  condition  a>x  +  coy  =  <u-.  For  sake  of  argu¬ 
ment,  we  assume  this  resonance  condition  to  be  exactly  satisfied.  We  consider  heating  of 
the  x  motion  assuming  both  the  y  and  z  modes  are  excited.  From  Eq.  (26),  we  find 
\x(t)\  =  \qtAyAz[d2Ex/dy  dz](JI=z=ty/{2ma>x)\  where  A v  and  A,  are  the  amplitudes  of  motion 
in  the  y  and  z  directions.  Here,  we  neglect  the  fact  that  energy  is  exchanged  between  the 
three  modes;  for  simplicity  we  assume  the  amplitudes  of  the  y  and  z  motion  remain  fixed. 
In  this  approximation,  if  Ay  =  AT  =  (;,  the  time  it  takes  to  excite  the  x  motion  to  the  same 
amplitude  is  given  by  t  =  \2mwx/ (q£\d2Ex/ dy  dz\^y=,  =  0)\.  If  £  =  10  nm  (corresponding  to 

(ny)  ~  («-)  ~  1  for  the  conditions  of  the  single  9Be+  ion  NIST  experiments,  the  field 
gradient  required  to  drive  the  x  motion  to  an  amplitude  of  10  nm  ((nx)  ~  1)  in  the 
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observed  time  of  1  ms  is  approximately  d2Ex/dydz=  1000  V/mm3.  It  is  highly  unlikely 
the  gradient  was  this  large  for  the  NIST  experiments,  and,  furthermore,  the  resonance 
condition  was  only  approximately  satisfied.  Moreover,  this  source  of  heating  was  easily 
tested  by  varying  the  initial  values  of  Ay  and  A,  (by  varying  the  Doppler-cooling  mini¬ 
mum  temperature  through  laser  detuning)  and  studying  the  heating  rate  of  the  x  motion 
which  had  previously  been  cooled  to  the  zero  point  of  motion.  No  dependence  on  the 
initial  values  of  Ay  and  Az  was  found.  In  any  case,  if  all  modes  of  motion  are  initially 
cooled  to  the  zero-point  state  this  source  of  heating  vanishes  because  the  assumed  cou¬ 
pling  only  provides  an  exchange  of  energy  between  modes.  This  places  a  premium  on 
cooling  all  modes  to  as  low  an  energy  as  possible.  Finally,  it  appears  that  this  single¬ 
ion  example  gives  a  worst  case  scenario  since,  for  large  numbers  of  ions,  the  force  on 
the  generalized  coordinates  (right  hand  side  of  Eq.  (23))  requires  a  high-order  field  gra¬ 
dient  to  be  nonzero.  These  gradients  are  highly  suppressed  in  the  typical  case  where  ion-ion 
separation  is  much  smaller  than  the  distance  between  the  ions  and  the  trap  electrodes. 

3.  Mode  cross  coupling  induced  by  logic  operations 

In  the  preceding  sections,  we  have  assumed  that  when  transitions  are  driven  between 
\nk)  |  j}7-  and  \n'k)  ||);  involving  a  single  mode  of  motion  k,  the  other  3/V  —  I  spectator 
modes  of  motion  are  not  affected  because  coupling  to  them  is  nonresonant.  However,  when 
the  sum  or  difference  frequency  of  two  or  more  spectator  modes  is  near  the  frequency  of 
the  desired  mode-k  transition  (wk  \n'k  —  nk\),  higher  order  couplings  can  entangle  the 
\nk)  |  j);  and  | n'k)  ||)  -  states  with  the  spectator  mode  states. 

Equation  (19)  describes  the  general  interaction  Hamiltonian  between  the  internal  levels  of 
ion  j  and  all  3N  modes  of  motion.  By  expanding  the  exponential  in  Eq.  (19)  to  all  orders, 
we  find 


H\  =  HQ' 


Hfli)  '  '  jQn-ddaA  J(d i- 

b,W  j 


h.c. 


(27) 


where  Qj  =  Q1  exp  [—  i  H/(?/,)2].  This  equation  describes  the  processes  of  each  mode  /  gain¬ 
ing  or  losing  ( bj  —  d/)  vibrational  quanta,  accompanied  by  the  raising  or  lowering  of  the 
internal  electronic  levels  of  ion  j.  We  must  account  for  all  terms  in  Eq.  (27)  which  do  not 
vary  rapidly  in  time,  or  terms  in  which  the  resonance  condition  is  nearly  met: 
S/(bi  —di)(o i  ~  (5  =  0)k{n'k  — nk ).  Although  detailed  treatment  of  this  problem  is  beyond 
our  intent,  a  couple  of  comments  may  be  made. 

In  general,  we  must  account  for  all  the  terms  in  Eq.  (27)  which  cause  significant  errors 
in  the  overall  computation  we  are  trying  to  carry  out.  For  two  or  more  trapped  ions,  some 
combination  of  modes  will  nearly  always  satisfy  the  resonance  condition.  However,  this 
may  occur  only  for  high  orders  of  b/  and  d /,  and  if  the  Lamb-Dicke  criterion  is  met.  the 
contributions  are  vanishingly  small.  The  terms  that  will  cause  problems  are  the  ones  that 
are  close  to  satisfying  the  resonance  condition  and  are  relatively  low  order  in  bj  and  d/.  If 
the  Lamb-Dicke  criterion  is  satisfied,  it  will  always  be  possible  to  avoid  these  spurious 
couplings,  but  it  may  be  at  the  expense  of  making  the  Rabi  frequency  so  small  (in  order  to 
avoid  coupling  to  relatively  nearby  off-resonant  terms)  that  the  operations  will  become  too 
slow. 

To  understand  this  problem  in  the  context  of  a  simple  example,  we  assume  that  a  cross¬ 
mode  coupling  of  this  type  occurs  when  two  modes,  p  and  q,  have  frequencies  which 
satisfy  the  condition  npa>p  —  nqcoq  ~  0,  a>k,  or  —  mk  corresponding  to  possible  extraneous 
mode  coupling  on  the  carrier,  first  blue  sideband,  or  first  red  sideband  of  the  logic  opera¬ 
tions  (assumed  to  utilize  mode  k ).  This  additional  resonance  condition  yields,  to  lowest 
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order  in  the  Lamb-Dicke  parameters,  the  effective  Hamiltonian 

H{  =  fi&jSi  1 1  +  irfk(ak  +  a\  ei0>tt) 

+  (»7 ,4TP  e,H 

np\nq\ 

where  the  resonance  conditions  are  given  by  d  =  mL  —  at0  =  0,  +o>k,  or  —<»k.  A  specific 
example  is  relevant  to  the  NIST  single  yBe  ion  experiments.  Here,  mode  k  was  the  x 
oscillation,  and  modes  p  and  q  are  identified  with  the  z  and  y  oscillations  of  the  single  ion  in 
the  trap.  In  this  experiment  i]x  =  dk  •  jf.ro,  Vy  =  dk  •  j(vo,  rjz=  dk  •  zzo ,  and  ojx  ~  to.  —  a>y. 
(The  frequency  relationship  a>x  =  o>-  —  a>y  is  a  consequence  of  Maxwell’s  equations  for  a 
quadrupole  rf  trap  in  the  absence  of  static  potentials  applied  to  the  electrodes  [43])  We 
assume  that  the  desired  transition  is  the  first  blue  sideband  of  mode  x(d  =  a)x).  In  this 
case,  the  resonant  part  of  Eq.  (28)  becomes 

H\  =  hQ'j[S]+{iqxa\  -  r]zr]ya[ay  +  0(t/3)}  +  h.c.] .  (29) 

The  term  proportional  to  a\  is  the  desired  anti-Jaynes-Cummings  operator,  and  the  term 
proportional  to  a[ay  can  entangle  the  internal  state  with  the  other  spectator  modes  (z  and  y), 
leading  to  errors. 

For  logic  operations  on  a  string  of  ions  in  a  linear  trap,  we  will  assume  that  all  other 
mode  frequencies  are  higher.  With  the  use  of  stimulated-Raman  transitions,  we  can  make 
<3Ar  ||  £  and  restrict  our  attention  to  spectator  modes  along  the  z  axis.  Nevertheless,  as  N 
becomes  large,  nearby  resonances  of  the  type  shown  in  Eq.  (28)  will  become  harder  to 
avoid.  These  coupling  terms  always  scale  as  products  of  Lamb-Dicke  parameters.  Thus  if 
the  spectator  mode  Lamb-Dicke  parameters  are  small  enough,  or  if  at  least  one  Lamb-Dicke 
parameter  is  approximately  zero,  the  higher  order  unwanted  resonances  may  be  sufficiently 
suppressed.  Furthermore,  if  the  spectator  modes  are  cooled  to  near  the  zero-point  energy 
((«)  1),  then  any  couplings  in  Eq.  (28)  with  powers  of  the  annihilation  operator  aq  will 

be  absent  most  of  the  time.  Hence,  in  large  registers,  it  will  probably  be  important  to  cool 
all  modes  to  near  the  zero-point  energy. 


-K8t-<pj)+  hc. 


(28) 


V.  Decoherence 

A.  Internal  state  decoherence  from  spontaneous  emission 

The  internal  atomic  states  of  trapped  ions,  which  store  quantum  bits  of  information,  must 
be  protected  from  spontaneous  emission,  at  least  for  the  duration  of  the  computation.  This 
excludes  the  possibility  of  “error  correction,’’  [48]  which  can  tolerate  a  certain  level  of 
errors  due  to  spontaneous  emission.  For  qubit  levels  coupled  by  single  photon  optical  tran¬ 
sitions,  this  may  be  accomplished  by  employing  long-lived  energy  levels  which  do  not  have 
an  allowed  electric  dipole  coupling,  such  as  metastable  electronic  levels  with  a  quadrupole 
or  intercombination  coupling  to  the  ground  state.  Although  the  interaction  of  these  states 
with  the  vacuum  (causing  spontaneous  emission)  is  reduced,  their  interaction  with  an  exter¬ 
nal  field  for  use  in  quantum  logic  operations  is  also  reduced.  This  results  in  a  fundamental 
limit  on  the  accuracy  of  each  operation  by  roughly  the  ratio  of  the  spontaneous  emission 
rate  to  the  Rabi  frequency  £  =  T/Q.  In  the  case  of  optical  transitions,  Q  cannot  be  in¬ 
creased  indefinitely,  since  at  optical  intensities  beyond  about  1014  W/cnr,  the  atom  is 
quickly  photoionized.  This  amounts  to  inaccuracies  due  to  spontaneous  emission  on  the 
order  of  £  =  10  6  — 10  7  [49].  Even  this  limit  may  be  too  optimistic,  as  the  two-level 
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approximation  breaks  down  before  photoionization  occurs,  and  the  coupling  to  other  elec¬ 
tronic  levels  must  also  be  considered  [49].  This  results  in  inaccuracies  due  to  spontaneous 
emission  on  the  order  of  £  =  10  s  — 10~6,  depending  on  the  particular  ion  used. 

In  the  case  of  two-photon  stimulated-Raman  transitions  between  stable  electronic  ground 
states,  the  ratio  of  spontaneous  emission  rate  to  Rabi  frequency  is  approximately 
£sr  =  Y  / (g2 / 4l),  where  yse  w  Fg1  /A1  is  the  off-resonant  spontaneous  emission  rate,  g  is 
the  resonant  single-photon  Rabi  frequency  of  each  laser  beam,  and  A  is  the  detuning  of  the 
Raman  beams  from  the  excited  state.  This  results  in  an  inaccuracy  r /A  due  to  spontaneous 
emission,  which  is  independent  of  optical  intensity.  Since  Raman  transitions  between  S 
electronic  ground  states  are  effective  only  when  the  detuning  A  is  not  much  greater  than 
the  fine  structure  splitting  of  the  atom  [50];  this  results  in  an  inaccuracy  £SR  due  to  sponta¬ 
neous  emission  in  range  from  about  10  4  (9Be+)  to  10  7  (199Hg+),  depending  on  the  parti¬ 
cular  ion  used.  Spontaneous  emission  from  spectator  electronic  levels  should  not  signifi¬ 
cantly  affect  this  limit,  provided  that  their  splitting  from  the  virtual  excited  state 
significantly  exceeds  A  and  that  the  single  photon  resonant  Rabi  frequencies  coupling  the 
ground  states  to  the  spectator  levels  are  not  much  bigger  than  g  [49].  These  appear  to  be 
reasonable  assumptions  for  most  candidate  ions. 

The  decohering  effects  of  spontaneous  emission  can  be  overcome  by  error  correction 
schemes.  Error  correction  is  complicated  by  the  fact  that  when  spontaneous  emission  oc¬ 
curs,  the  atoms  may  decay  to  states  outside  the  original  set  of  computational  basis  states. 
However,  this  situation  can,  in  principle,  be  detected  by  optically  pumping  the  ions  back  to 
the  computational  basis  and  applying  the  error  correction  schemes  [49,  51]. 

Spontaneous  emission  decoherence  could,  in  principle,  be  nearly  eliminated  by  driving 
single-photon  transitions  between  ground-state-hyperfine  or  Zeeman  levels  with  rf  or  micro- 
wave  radiation  since  spontaneous  emission  from  these  levels  is  negligible.  This  may  be 
accomplished  by  coupling  the  internal  and  motional  states  with  inhomogeneous  magnetic 
fields  [52].  The  speed  of  sideband  operations  is  limited  by  the  size  of  the  field  inhomogene¬ 
ity  one  can  obtain.  For  example,  consider  an  ion  moving  along  the  ’-axis  with  a  magnetic 
field  gradient  dBx/dz  applied  on  top  of  a  uniform  magnetic  field  B  along  the  direction  of 
motion.  In  the  ion’s  reference  frame,  there  is  an  rf  magnetic  field  which  can  induce  transi¬ 
tions  between  the  internal  states.  The  interaction  Hamiltonian  is  then 

n  =  -fi  B  =  —fiBzSz  -  /'(dftr/fcho  +  s_^  (flt  +  a) ;  (30) 

where  u  «  iiB  is  the  magnitude  of  the  ion  magnetic  moment,  typically  near  one  Bohr  mag¬ 
neton,  ’o  is  the  spatial  spread  of  the  zero-point  ion  wavepacket,  S±  are  the  spin  raising  and 
lowering  operators  from  Sec.  II.C  and  a  and  «  are  the  usual  motional  raising  and  lowering 
operators.  The  coupling  (last  term  in  Eq.  (30))  is  analogous  to  Eq.  (3),  allowing  the  con¬ 
struction  of  logic  gates  as  described  above.  Unfortunately,  this  method  requires  very  high 
magnetic  field  gradients  dBx/dz.  For  instance,  in  order  to  realize  an  effective  Rabi  fre¬ 
quency  of  0.1  MHz  (£2eff  ~  (A B(dBx / dz)  Zo / 4) ,  a  magnetic  field  gradient  of  around  10  T/cm 
would  be  necessary,  where  we  have  assumed  that  the  motional  frequency  is  1  MHz. 
Although  these  gradients  aren’t  unreasonable  in  the  laboratory,  the  residual  fields  that  ac¬ 
company  the  gradient  will  dramatically  shift  the  internal  electronic  levels  of  the  trapped 
ions.  Moreover,  it  would  be  difficult  to  address  selected  ions  because  of  the  long  wave¬ 
length  of  the  radiation  relative  to  typical  ion  spacings. 

B.  Motional  decoherence 

Although  the  internal  states  of  ions  may  be  well  isolated  from  environmental  influences,  the 
motional  states  of  a  collection  of  trapped  ions  are  expected  to  be  more  susceptible  to  deco- 
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herence.  The  Cirac/Zoller  scheme  for  quantum  logic  [9]  transfers  quantum  information  tran¬ 
siently  through  the  motional  states,  so  any  decoherence  of  the  motion  will  compromise 
operation  of  any  multi  bit  quantum  gate  and  must  be  minimized.  In  the  NIST  experiments 
with  a  single  ion,  the  motional  state  of  all  three  modes  was  observed  to  absorb  energy  from 
the  environment  at  the  rate  of  about  1  phonon  (energy  ha)z,  where  (o./2ji  ss  10  MHz)  per 
millisecond  from  the  ground  state  [10,  17].  Although  this  as  of  yet  unexplained  heating  rate 
is  not  believed  to  be  fundamental,  it  is  of  considerable  interest  to  characterize  the  possible 
sources.  There  are  several  potential  mechanisms  for  motional  decoherence;  for  instance  “rf 
heating”  from  rf  fields  in  the  trap  [53],  collisions  with  background  atoms,  fluctuating  patch 
and  contact  potentials  on  the  trap  electrodes,  electron  field  emission  from  the  electrodes, 
and  Ohmic  loss  of  induced  image  currents  in  the  trap  electrodes  [54,  55].  These  and  other 
effects,  considered  in  Ref.  [11],  are  not  expected  to  be  major  obstacles  to  motional  coher¬ 
ence  in  a  linear  ion  trap. 

By  virtue  of  the  electric  charge  of  trapped  ions,  it  is  natural  to  consider  the  coupling  of 
spurious  electric  fields  to  the  motion  of  the  ion  string.  In  this  section  we  relate  the  size  of 
electric  field  noise  to  the  expected  transition  rate  from  the  quantum  ground  state  of  motion 
in  the  ion  trap.  Such  noisy  electric  fields  might  be  generated  from  unstable  trap  parameters, 
external  radiation,  thermal  (Johnson)  noise  from  resistive  losses  in  the  trap  electrodes,  or 
patch  potentials  from  nearby  surfaces. 

We  first  consider  the  axial  COM  harmonic  motion  of  a  string  of  N  ions  at  frequency  a>z. 
An  external  uniform  electric  field  will  shift  the  position  of  the  trap  center;  an  external 
electric  field  gradient  will  change  the  effective  spring  constant  of  the  trap.  Following  the 
approach  of  Savard,  O’Hara  and  Thomas  [56],  who  considered  heating  mechanisms  in 
the  context  of  neutral  atom  dipole  traps,  we  can  similarly  calculate  the  rate  C 0  at  which  the 
ground  state  of  ion  motion  is  vacated  due  to  noisy  electric  fields.  We  find 


/-'axial 

1  0 


Nq2 

4mficoz 


Se{ojz) 


H 

4Nma)z 


Se'(2ojz)\  , 


(31) 


where  Se(cl>)  is  the  electric  field  noise  spectral  density  in  (V/cm)2/Hz  and  SE'(a>)  is  the 
electric  field  gradient  noise  spectral  density  in  (V/cm2)2/Hz.  The  noise  densities  are  defined 
so  that  the  mean-squared  electric  field  and  field  gradient  are  (2jt)_1  da)  Se(oj)  and 
(2jt)_1  J0°°  dajSE'(a >),  respectively  [56]. 

Note  that  contributions  of  electric  field  noise  to  motional  heating  is  concentrated  near 
resonance  (co.)  in  Eq.  (31),  as  expected  for  a  resonantly  driven  harmonic  oscillator.  Contri¬ 
butions  of  electric  field  gradient  noise  to  heating  is  concentrated  near  twice  the  resonant 
frequency  (2a>z),  as  expected  for  a  parametrically  driven  harmonic  oscillator.  For  modes  of 
motion  with  secular  frequency  a>x  derived  from  ponderomotive  rf  forces  (i.e.,  the  radial 
modes  of  a  linear  trap),  there  are  additional  contributions  to  heating  from  noisy  fields  at 
frequencies  Qj  ±  wx  and  gradients  at  frequencies  Qj  ±  2ojx,  to  lowest  order  in  the  pseudo¬ 
potential  approximation  \(wx/Qt)  1]  [21].  We  find 


Se((jdx)  H - \  Se(£2t  3^  tux)  3 - [  S f  ( 2ojx )  -I - A  SE/(^2E  i  2cuY)  |  . 

2  Q].  4Nmcoxy  V  x>  f22  ) 

(32) 

Given  a  source  of  electric  field  noise  Se(<w),  we  may  expect  the  corresponding  gradi¬ 
ent  noise  spectrum  S^cu)1'2  to  be  no  larger  than  Se{(d)1  2 /d,  where  d  is  the  characteristic 
ion  trap  electrode  size.  From  the  above  equations,  we  find  that  this  implies  that  the  heating 
rate  from  gradient  noise  is  expected  to  be  a  factor  of  «  ( zo/2d )2  smaller  than  the  heating 
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rate  from  electric  field  noise,  where  we  recall  that  zo  =  (. ti/2Nmcoz )*'2  is  the  size  of  the 
ground  state  wavepacket  of  the  mode  of  interest.  This  ratio  is  typically  very  small  in  ion 
traps  (in  the  NIST  9Be~  experiments,  zo  ~  7  nm,  d  «  250  |im,  and  (z,o/2d)~  k2x  10~10). 
Noise  in  the  electric  field  gradient  is  therefore  expected  to  be  much  less  of  a  concern  than 
noise  in  the  electric  field,  unless  the  electric  potentials  on  the  trap  electrodes  possess  an 
unusually  high  degree  of  symmetry.  It  also  follows  that  the  COM  modes  should  be  most 
susceptible  to  heating,  as  the  non-COM  modes  do  not  respond  to  uniform  electric  fields. 
For  instance,  in  the  case  of  two  ions  confined  along  the  z  axis,  the  COM  mode  along  z  is 
sensitive  to  resonant  uniform  electric  fields.  The  “stretch”  mode  along  z  (where  the  two 
ions’  motions  are  out  of  phase)  on  the  other  hand,  is  only  sensitive  to  resonant  electric  field 
gradients.  Therefore,  if  quantum  logic  operations  are  limited  by  motional  heating  due  to 
uniform  electric  fields,  higher-order  (non-COM)  motional  modes  should  be  used  for  quan¬ 
tum  logic.  COM  heating  will  still  indirectly  affect  the  fidelity  of  operations  on  other  modes, 
but  this  effect  should  be  of  higher  order,  as  discussed  in  section  IV.B.l. 

We  examine  two  expected  sources  of  electric  field  noise,  and  estimate  their  effects  on  the 
heating  rate  of  COM  modes.  These  results  are  compared  to  the  NIST  experiment  on  a 
single  9Be+  ion  [17],  where  we  observed  a  heating  rate  To  «  103/sec,  corresponding  to  an 
electric  field  noise  of  Se{u>x)X^2  ~  4  •  10~8  V/cm/Hz1/2  or  electric  field  gradient  noise  of 
SE'(2cox)1/2  «  0.07  V/cnr/Hz1/2. 


1.  Thermal  or  blackbody  noise 

Thermal  fluctuations  from  lossy  elements  of  the  trap  electrode  structure  may  lead  to  electric 
field  noise  at  the  ions.  This  is  essentially  the  effect  of  blackbody  radiation  on  the  ion 
motion,  altered  by  the  “cavity”  formed  from  the  trap  electrodes  [11,  54,  57].  Lossy  or 
resistive  elements  of  the  trap  electrodes  will  only  give  rise  to  a  noisy  electric  field  only  if 
the  effective  resistive  current  path  is  asymmetrically  oriented  with  respect  to  the  ions;  other¬ 
wise  only  neglegible  field  gradients  will  contribute,  as  discussed  above.  The  voltage  noise 
spectral  density  from  a  resistor  R  at  temperature  T  is  AkTR.  For  an  electrode  structure  of 
characteristic  size  d,  we  thus  expect  an  electric  field  noise  density  of  Se((o)  «  AkTR{co)/d 2, 
where  R(a>)  describes  the  resistance  seen  by  the  ion  from  the  surrounding  electrode  struc¬ 
ture,  or,  equivalently,  the  resistance  through  which  currents  induced  by  ion  motion  flow. 
Neglecting  the  gradient  terms  in  Eqs.  (31)  and  (32),  this  results  in  axial  and  radial  mode 
heating  rates  of  [7,  11] 
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We  now  use  this  last  equation  to  estimate  the  heating  rate  which  might  have  been  ex¬ 
pected  in  the  NIST  single  ion  9Be+  experiments  [17],  where  all  modes  were  confined  from 
pondero motive  forces.  We  consider  two  current  paths  in  the  electrode  structure  which 
would  likely  have  provided  the  largest  resistances.  The  first  is  a  direct  path  between  the 
endcap  electrodes,  which  were  positioned  around  the  ring  electrode  [11,  43].  Since  this 
path  length  is  only  about  / «  0.5  cm  and  much  smaller  than  the  wavelengths  associated 
with  the  frequencies  of  interest  mx/2n  «  10  MHz  and  (Qj  ±  mx)/2jz  «230  MHz,  we  can 
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treat  this  path  as  a  lumped  resistor  with  resistance  R(a>)  w  gl/A(co).  Here,  g  is  the  resistiv¬ 
ity  of  the  electrode  material  and  A(co)  is  the  effective  cross-sectional  area  of  the  resistive 
path,  proportional  to  the  skin  depth  at  frequency  to.  We  conservatively  estimate 
R(a>x)  «  0.04  Q  and  R(Q j  ±  wx)  «  0.20  Q  in  the  experiment  [11,  43]  and  calculate  a  heat¬ 
ing  rate  of  T0  «  0.7/s,  dominated  by  the  first  term  of  Eq.  (34)  and  much  smaller  than 
observed. 

The  rf  and  static  electrodes  in  an  rf  trap  are  typically  joined  through  an  rf  step-up  trans¬ 
former.  This  leads  us  to  consider  a  resistive  path  from  the  rf  to  the  static  electrode,  as  the 
resistance  between  these  electrodes  can  be  high  at  frequencies  near  the  resonant  rf  drive 
frequency  QT.  In  a  linear  trap,  these  fields  should  not  have  a  component  along  the  axial 
modes  of  motion;  moreover,  these  fields  will  vanish  at  trap  center  and  only  provide  a  gradi¬ 
ent,  unless  there  is  an  extreme  asymmetry  in  the  electrode  structure.  We  thus  consider  noise 
fields  resulting  from  a  geometrical  asymmetry  ft  between  the  ring  and  endcap  electrodes 
(0  <  ft  <  1,  where  ft  is  proportional  to  the  electric  field  at  the  ion  if  a  potential  is  applied 
between  the  ring  and  endcaps;  ft  =  0  lor  perfect  symmetry).  In  the  NIST  9Be+  experiments 
[10,  17],  a  quarter-wave  transmision  line  acts  as  the  step-up  transformer  [43],  thus  the 
impedance  between  the  rf  and  static  trap  electrodes  as  a  function  of  frequency  a>  is 


Z.rapM  =  Z(i  tanh 


jt 
4 Q 


JT 

2. Q~r 


(35) 


where  Zq  is  the  characteristic  line  impedance,  and  Q  is  the  loaded  quality  factor  of  the 
transmission  line  of  length  jic/2Qt  at  resonance  (a>  =  Q T).  The  resistance  is  the  real 
part  of  Eq.  (35),  and  at  the  frequencies  of  interest  which  might  couple  to  the  ion  mo¬ 
tion,  a>x  and  Qj  ±  a)x,  we  have  R(a>x)  «  (jtZq/AQ)  (cJx/Qt)1^  m  0.03  Q  and 
R(Qt  ±  coft)  «  (Zq/jtQ)  (Qt/o>x)2  s=s  34  £2.  In  these  expressions,  Zq  «  100  Q  and  Q  «  500 
are  the  characteristic  impedance  and  quality  factor  of  the  transformer,  and  (dx/2jt  and 
Qt/2 jt  are  taken  to  be  10  MHz  and  230  MHz  respectively  [43].  Note  the  estimated  resis¬ 
tance  R(a>x)  is  nearly  equivalent  to  the  value  in  the  skin-depth  model  above,  since  this  path 
length  is  still  only  a  small  fraction  of  the  wavelength  associated  with  cox.  We  find  that  the 
two  terms  in  Eq.  (34)  contribute  roughly  equal  amounts  to  a  total  heating  rate  of  approxi¬ 
mately  To  «  0.8/T/s.  Again,  this  rate  is  much  smaller  than  the  observed  rate,  even  for 
large  asymmetries. 


2.  Noise  on  trap  voltages 

Fluctuations  in  the  trap  voltages  (t/o  and  V0  in  Eq.  (1))  will  nominally  only  give  rise  to  a 
noisy  field  gradient.  However,  in  some  cases,  asymmetric  fluctuations  in  trap  parameters 
can  give  rise  to  a  noisy  electric  field.  For  instance,  noise  on  the  static  potentials  of  the  end 
segments  of  a  linear  trap  (Fig.  2)  may  not  be  common  mode  if  the  potentials  are  provided 
from  uncorrelated  power  supplies,  or  if  the  leads  connecting  to  the  end  segments  experience 
different  amounts  of  noise  from  pickup.  In  these  cases,  we  can  relate  the  electric  field  noise 
density  to  a  differential  potential  noise  density  defined  by  Sm0(a>)  =  d2SE(a>),  where  dUo 
is  the  difference  in  potentials  between  the  rod  end  segments.  From  Eq.  (31),  this  would 
give  rise  to  a  heating  rate  of 


Nq 2 

4  mfia):d2 


SdUoiWz)  ■ 


(36) 


For  example,  a  typical  power  supply  might  have  rms  noise  of  0.1  mV  uniformly  distributed 
across  a  20  MHz  bandwidth.  If  the  supply  is  filtered  so  that  signal  amplitude  at  frequency 
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cux  is  attenuated  by  a  factor  F  <  1,  we  find  [S,5(/0(2a>,)]1'/”  «  20 F  (nV/Hz1/2),  leading  to  a 
heating  rate  of  the  axial  COM  mode  of  Co  ss  4  •  106/V  •  F2 / s.  Here  we  have  assumed  trap 
parameters  (oz/2jt  ~  1  MHz  and  d  «  0.025  cm. 

If  there  exist  asymmetric  static  patch  or  contact  potentials  on  the  electrodes,  this  could 
lead  to  a  static  electric  field  fsSLa,;c  which  would  push  the  ions  away  from  at  the  geometrical 
trap  center  defined  byx  =  y  =  z  =  0  in  Eq.  (1).  In  this  case,  fluctuations  in  the  potentials 
U0  and  Vo  will  be  converted  into  a  noisy  electric  field  at  the  ions.  If  Est atiC  is  in  the  axial 
direction,  then  a  noise  density  of  the  static  voltage  Su0(co)  would  be  equivalent  to  a  field 
noise  density  of  Se(o)  =  (£’static/(Co))2  Su0(a>).  If  we  assume  Zsstatjc  ~  (p/d,  where  cp  is  an 
effective  axial  potential  difference  across  the  electrodes  (due  to  patch  potentials  for  in¬ 
stance),  then  we  find  from  Eq.  (31) 


Nq2 

4mHa>zd 2 


(37) 


Similar  to  above,  this  leads  to  a  heating  rate  of  .To  ~  4  •  106(V  •  F2(cp/(Uo))2 /s  assum¬ 
ing  the  above  values  of  the  trap  parameters.  Controlling  these  potentially  troublesome 
heating  rates  clearly  points  to  the  importance  of  heavily  filtering  the  static  electrodes  at 
frequency  wz.  The  filtering  is  best  accomplished  as  close  as  possible  to  the  trap  electro¬ 
des,  and  for  reasonable  filter  factors  of  F  <  10  4,  these  heating  sources  should  be  neg¬ 
ligible. 

For  fluctuations  in  the  rf  potential  Vo,  we  can  derive  a  similar  heating  rate  of  motional 
modes  confined  by  ponderomotive  rf  forces,  which  would  be  relevant  to  the  NIST  experi¬ 
ments.  For  simplicity,  we  consider  the  COM  radial  mode  at  frequency  cox  and  assume  the 
confinement  is  dominated  by  the  ponderomotive  force  (cop  term  in  Eq.  (2)).  We  consider  the 
case  of  a  patch  field  £static  =  <p/R-  The  equivalent  electric  field  noise  density  is  then 
Se(o>)  =  (2Estatic/(Vo))2  Sv0 (to),  which  from  Eq.  (32)  results  in  a  heating  rate 


/-radial 
1  0 
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where  Sy0(a>)  refers  to  the  rf  potential  noise  density  at  the  trap  rf  electrode.  To  estimate  the 
effect  of  rf  amplitude  noise  at  the  input,  we  again  assume  that  the  rf  input  lead  is  connected 
to  a  step-up  transformer.  For  a  transformer  of  characteristic  impedance  Zq  and  quality  factor 
Q  »  Qj/2mx  >>  1,  input  signals  of  frequency  <nx  are  essentially  filtered  out,  and  signals  at 
Qt  ±  ojx  are  multiplied  by  the  factor  [Zq/ (2tcRsQ)\  12  Qt/ojx,  where  Rs  is  the  source  impe¬ 
dance.  We  thus  neglect  the  first  term  of  Eq.  (38)  and  replace  Sv0{&t  ±  (ox)  by 
[Zo/{2jtRsQ)\  ( Qt/wx)2SVo ,  input  (^r  ±  (ox)  in  the  second  term.  If  the  rf  source  feeding  the 
input  has  an  effective  noise  figure  of  NF  «  10  dB  above  the  Johnson  noise  of  a  Rs  =  50  £2 
source  impedance,  we  have  SV„, input  ~  lO(4kTRs)  w  1CV17  V2/Hz.  Assuming 
cox/2 jt «  10  MHz,  we  find  (Vo)  w  750  V  and  a  heating  rate  of  To  «  10(<)9/(Vo))2/s.  Since 
(p  (Vo)  for  typical  patch  potentials  of  <  1  V,  this  source  of  heating  is  again  much  smal¬ 
ler  than  observed.  Moreover,  in  the  NIST  experiments,  cp  was  varied  over  a  wide  range 
without  an  observed  dependence  on  heating  rate. 

Although  we  have  tried  to  characterize  the  more  obvious  sources  of  electric  field 
noise,  the  actual  electric  field  noise  spectral  density  in  an  ion  trap  may  have  a  compli¬ 
cated  structure  and  may  be  difficult  to  characterize.  The  motional  modes  of  a  string  of 
ions  are  most  susceptible  to  electric  field  noise  near  the  motional  frequency  of  interest, 
as  expected.  If  spurious  resonances  should  occur  at  these  frequencies  due  to  background 
electric  fields,  it  might  be  possible  to  avoid  them  by  simply  changing  the  trap  param- 
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eters.  As  stated  above,  heavy  filtering  of  the  electrodes  should  provide  adequate  noise 
suppression  at  these  frequencies  and  the  COM  motion  of  ions  trapped  in  electrode 
structures  with  a  high  degree  of  symmetry  will  be  less  susceptible  to  electric  field 
noise. 


C.  Induced  decoherence  from  applied  field  amplitude  noise 


In  the  ion  trap  quantum  computer,  all  operations  can  be  traced  back  to  Eq.  (5),  the  coupling 
of  states  | n)  |j)  to  | n')  ||).  The  fidelity  of  these  operations  will  depend,  in  part,  on  how 
accurately  the  coupling  strength  Q1  and  the  application  time  can  be  set.  For  instance,  in  the 
rotation  gate  of  Eq.  (7),  the  rotation  angle  0  =  2f2/f,  and  noise  in  these  parameters  will  lead 
to  evolution  to  an  undesired  quantum  state.  The  effects  of  (Gaussian)  noise  on  laser  inten¬ 
sity  have  been  treated  by  Schneider  and  Milburn  [58].  These  effects  show  up  in  a  well- 
characterized  way  for  transitions  involving  Fock  states.  Here  we  briefly  investigate  the  size 
of  laser  power  fluctuations,  as  the  Rabi  frequency  is  proportional  to  the  laser  power  in  the 
Raman  configuration.  We  do  not  examine  how  such  errors  might  propagate  in  an  extended 
gate  structure  [9,  11,  59]. 

Fluctuations  in  the  laser  intensity  at  the  site  of  a  given  ion  can  arise  from  both 
fluctuations  in  the  relative  position  of  the  beam  with  respect  to  the  ion  and  fluctuations 
in  laser  power.  Laser/ion  position  stability  is  particularly  important  since  the  Cirac/Zoller 
scheme  of  quantum  logic  assumes  that  ions  in  an  array  be  selectively  addressed,  there¬ 
by  requiring  a  high  degree  of  control  of  the  laser  beam  spatial  profile  (Sec.  IV.A).  Of 
course,  the  simplest  method  for  minimizing  position  fluctuations  is  to  employ  mechani¬ 
cally  stiff  mounts  for  the  optics  and  ion  trap  electrodes,  and  have  the  laser  source  as 
close  as  possible  to  the  ions.  A  quadrant  detector  indexed  to  the  trap  electrodes  and 
placed  near  the  ion  may  also  be  used  to  actively  stabilize  the  beam  position  by  feeding 
back  to  a  galvanometer  or  acousto-optic  modulator.  If  optical  fibers  are  used  to  deliver 
laser  beams  to  the  ions,  position  fluctuations  between  the  fiber  and  the  ions  could  be 
made  small;  however,  we  must  also  consider  position  fluctuations  between  the  laser 
source  and  the  input  to  the  fiber.  If  the  position  tolerances  can  be  adequately  con¬ 
trolled,  the  dominant  source  of  intensity  noise  at  the  ion  would  likely  be  given  by 
fluctuations  in  optical  power  and  laser  mode.  Here,  we  estimate  limits  on  laser  ampli¬ 
tude  noise. 

If  we  assume  the  laser  fields  responsible  for  quantum  logic  operations  are  coherent 
states,  the  fundamental  noise  floor  is  photon  shot  noise.  For  a  laser  beam  of  average  power 
Pq  the  fractional  level  of  shot  noise  is 


dP  I  Tim 
Po  V  PoTop  ’ 


(39) 


where  co  is  the  (optical)  photon  frequency,  rop  is  the  time  the  radiation  is  applied,  and,  for 
simplicity,  we  assume  square  pulse  envelopes.  Almost  all  laser  sources  have  significant 
amplitude  noise  well  above  the  shot-noise  limit  in  the  10  Hz— 10  kHz  range  due  to  acoustic 
vibrations  which,  for  example,  affect  the  laser  cavity  resonators.  Much  of  this  noise  can  be 
removed  by  active  power  stabilization,  where  a  beamsplitter  directs  a  portion  of  the  laser 
power  to  a  photodetector,  and  an  error  signal  is  derived  and  fed  back  to  an  upstream  mod¬ 
ulator  or,  in  the  case  of  a  diode  laser,  directly  to  the  amplitude  of  the  laser  source  [60],  The 
limiting  noise  of  this  stabilization  scheme  is  degraded  slightly  by  the  imperfect  quantum 
efficiency  of  the  photodetector  as  well  as  the  beamsplitter.  If  the  beamsplitter  directs  a 
fraction,  of  the  input  optical  power  to  the  stabilizer  (which  then  gives  an  optical  power 
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Pu  ~  (1  —  e)  Pq  directed  to  the  ion),  and  the  quantum  efficiency  of  the  photo  detector  is 
j/det,  the  limit  of  fractional  power  noise  in  the  logic  pulse  is  (assuming  no  added  electronic 
noise  in  the  feedback  loop) 


dPu  >  I  km 

Pu  y  ^>HropJ/det£(l 


(40) 


For  a  laser  wavelength  of  313  nm,  and  assuming,  c  =  0.5,  and  rjdcl  =  0.5,  we  have 
6PJPU  >  2.3  ■  10_9(/>„rop)-1'2.  For  1  W  of  usable  laser  power  and  rop  =  1  [is,  this  corre¬ 
sponds  to  a  fractional  power  fluctuation  of  >  2.3  •  10  6. 

This  estimate  applies  only  to  the  laser  power  fluctuations  at  the  beamsplitter  and  assumes 
no  additional  noise  is  introduced  between  the  beamsplitter  and  the  photo  detector  or  the 
beam  splitter  and  the  ions.  Typically,  the  usable  part  of  the  laser  beam  must  be  directed 
further  through  optics,  the  air,  and  a  window  to  the  vacuum  envelope  enclosing  the  ion 
trap.  Fluctuating  etalon  effects  in  the  optics  and  air  currents  may  therefore  seriously  in¬ 
crease  the  power  fluctuations  beyond  Eq.  (39). 

Fluctuations  in  timing  errors  may  also  be  important.  Clearly,  fractional  fluctuations  in  the 
duration  of  laser  pulses  will  correspond  directly  to  the  same  fractional  fluctuations  in  the 
desired  value  of  the  gate  parameters  (e.g.,  rotation  angles  in  Eq.  (7)).  If  we  require  frac¬ 
tional  fluctuations  of  1 0  6  on  these  parameters,  then  we  require  timing  precision  of  1  ps  on 
a  1  s  pulse.  Similar  considerations  apply  to  the  stability  of  pulse  envelope  shapes. 

For  both  amplitude  and  timing  fluctuations,  it  may  be  possible  to  sample  a  portion  of  the 
beam  used  for  logic  and  apply  it  to  a  “check  bit”  ion.  The  response  of  this  ion  could  then 
be  used  to  monitor  and  control  the  amplitude  and  timing  of  the  pulses. 


VI.  Conclusion 

A  system  of  trapped  and  cooled  atomic  ions  is  one  of  the  few  viable  experimental  candi¬ 
dates  for  quantum  computation.  Internal  levels  of  the  ions  can  coherently  store  quantum 
bits  for  extremely  long  times.  Preserving  coherence  during  logic  gates  will  clearly  be  more 
difficult,  because  (1)  the  quantum  bits  are  distributed  through  a  collective  motional  degree 
of  freedom  in  the  trap  which  more  readily  couples  to  the  environment,  and  (2)  technical 
noise  on  the  logic  gates  degrades  gate  operation.  Nevertheless,  to  our  knowledge,  these 
difficulties  appear  to  be  technical,  not  fundamental.  The  maximum  number  of  gates  which 
can  be  applied  coherently,  and  the  maximum  number  of  ions  which  can  be  packed  into  a 
trap,  will  undoubtedly  be  determined  by  technical  limits.  Although  it  appears  that  the 
quantum  factoring  algorithm  would  be  extremely  difficult  to  implement,  we  conclude  that 
an  ion  trap  quantum  computer  of  very  modest  numbers  of  bits  and  gates  looks  quite  pro¬ 
mising. 

In  this  paper,  we  have  attemped  to  identify  some  of  the  more  important  experimental 
concerns  with  the  ion  trap  system,  and  expect  that  their  resolution  will  likely  be  determined 
in  the  laboratory.  We  summarize  with  a  few  general  observations  regarding  the  ion  trap 
quantum  computer. 

(1)  The  most  attractive  coupling  scheme  appears  to  involve  two-photon  stimulated  Ra¬ 
man  transitions  between  hyperfine  (or  Zeeman)  ground  states  separated  by  rf  or  microwave 
frequencies.  This  not  only  relaxes  the  laser  frequency  stability  requirements  when  compared 
to  single-photon  optical  transitions,  but  also  allows  more  geometrical  control  of  the  effec¬ 
tive  wavevector  dk  =  k\  —  ki.  which  determines  the  coupling  to  particular  motional  modes. 

(2)  Laser  cooling  of  all  modes  having  a  component  of  motion  along  dk  to  their  ground 
state  energy  will  minimize  cross  coupling  between  the  modes. 
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(3)  The  axial  COM  motional  mode  may  not  be  the  best  choice  for  a  “data  bus”  which 
couples  the  quantum  bits  between  ions.  First,  although  the  COM  mode  is  most  resolved 
from  the  other  modes,  it  has  the  lowest  motional  frequency,  which  may  ultimately  set  a 
limit  on  the  quantum  gate  speed.  Second,  the  non-COM  motional  modes  should  be  less 
sensitive  to  heating  from  external  electric  fields. 

(4)  For  a  given  number  of  ions,  the  optimum  value  of  the  mode  frequencies  will  likely 
involve  a  key  tradeoff.  Larger  frequencies  might  minimize  the  effect  of  some  sources  of 
motional  decoherence,  as  the  motional  energy  levels  would  be  further  separated.  This 
would  also  allow  gate  operations  (and  laser  cooling)  to  proceed  faster,  as  the  Rabi  fre¬ 
quency  can  be  no  larger  than  the  mode  frequency.  On  the  other  hand,  smaller  frequencies 
will  increase  the  spatial  separation  of  the  ions,  thus  easing  the  problem  of  optical  addres¬ 
sing  of  individual  ions  in  a  string. 

(5)  Internal  state  decoherence  will  likely  be  small,  and  it  does  not  appear  that  motional 
decoherence  will  be  a  fundamental  problem,  notwithstanding  the  heating  problem  observed 
in  the  NIST  experiments.  It  is  important  to  heavily  filter  the  ion  trap  static  electrodes  at  mz 
to  minimize  environmental  influences  on  the  ion  string  COM  motion. 
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